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Perturbative Renormalization 



1.1 Introduction 

Most of our present knowledge on the structure of renormalization of quantum field 
theories comes from perturbation theory (PT), a small coupling expansion around a 
free field theory. In this framework Feynman rules, derived formally in the continuum 
from the Gell-Mann-Low formula (fGcll-Man n and Low, 19 51) or from the path inte- 
gral approach, express amplitudes at a given order as sums of expressions associated 
with Feynman diagrams (see e.g. Nakanishi). Tree diagrams are associated with well 
defined amplitudes, but for diagrams involving loops and associated integration over 
internal momenta one soon encounters divergent expressions e.g. for massive go4> A 
theory in four euclidian dimensions the diagram in Fig. 11.11 is associated with the 
expression 

k {k 2 + m 2 ){[{k + p) 2 + m 2 ]) 

which is logarithmically divergent as the ultra-violet (UV) cutoff A —> oo. There 
are many ways to introduce an UV cutoff A - a specification defines a particular 
regularization. The process of renormalization is then a well-defined prescription how 
to map regularized expressions to amplitudes which are finite when A — > oo while 
maintaining desired properties which are summarized by the Osterwalder-Schrader 
axioms (jOsterwalder and Schrader 19731 119751) order by order in PT. 

which are equivalent to the Wightman axioms in Minkowski space (Wightman 1956 
|Streat cr and Wight man 1989fr 



k + p 




k 




Fig. 1.1 Simple 1-loop diagram 
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For lattice QCD renormalization is not needed if one is only interested in obtaining 
the low energy (LE) spectrum and scattering data. For such purposes we need knowl- 
edge of the phase diagram, the location of critical points where the continuum limit 
is reached, and the nature of the approach e.g. the question whether ratios of masses 
tend to their continuum limit as powers in the lattice spacing a: 

^ = ^§ ]+ C 12 (a mir , p>0(7). (1.2) 
m 2 (a) m 2 (0) 

Lattice artifacts will be the topic of chapters 4-6. But perturbative and non-perturbative 
renormalization will be needed for 1) computing matrix elements of composite op- 
erators (describing probes of other interactions, finite temperature transport coeffi- 
cients,....), 2) relating LE to high energy (HE) scales (e.g. computing running couplings 
and running masses), and 3) giving hints on the nature of lattice artifacts. 

In this chapter we shall mainly consider perturbative renormalization; in chap- 
ter 2 we shall discuss renormalization group equations which follow from multiplicative 
renormalization, and aspects of non-perturbative renormalization will be the subject 
of chapter 3. 



1.2 History and basic concepts 

In a series of papers Reisz (1988& 1989} has proven perturbative renormalizability of 
lattice Yang-Mills theory and QCD (for a large class of actions). His proof is based 
mainly on methods developed with continuum regularization, (with some important 
modifications that we will mention later), so we will start the discussion with these. 

Renormalization theory has a long and interesting history; here I reproduce a part 
of Wightman's delightful discussion ( [Wightman, 1975] ) since the reference is not always 
easily available. 

"One of the first steps in natural history is to establish a classifactory nomenclature. 
I will do this for perturbative renormalization theory, but in so doing, I want to tell 
stories with a moral for an earnest student. Renormalization theory has a history of 
egregious errors by distinguished savants. It has a justified reputation for perversity; 
a method that works up to 13th order in the perturbation series fails in the 14th 
order. Arguments that sound plausible often dissolve into mush when examined closely. 
The worst that can happen often happens. The prudent student would do well to 
distinguish sharply between what has been proved and what has been made 
plausible, and in general he should watch out! 

My first cautionary tale has to do with the early days of renormalization the- 
ory. When F.J. Dyson analyzed the renormalization theory of the S-matrix for quan- 
tum electrodynamics of spin one-half particles in his two great papers of 1948-9, 
{Dyson 19 49a, 1 949b]) he laid the foundations for most later work on the subject, but 
his treatment of one phenomenon, overlapping divergences was incomplete. Among 
the methods offered to clarify the situation, that of J. Ward {1951)) seemed outstand- 
ingly simple, so much so that it was adopted in Jauch and Rohrlich's standard text 
book {Jauch and Rohrl ich, 1955^. Several years later Mills and Yang noticed that un- 
less further refinements are introduced the method does not work for the photon self 
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Fig. 1.2 Troublesome 14'th oder QED diagram 

energy flWu, 19 62). The lowest order for which the trouble manifests itself is the four- 
teenth e.g. in the (7-loop) graph Fig. \1.2\ Mills and Yang repaired the method and 
sketched some of the steps in a proof that it would yield a finite renormalized am- 
plitude (Mills and Yang, 1966). An innocent reading of the textbook of Jauch and 
Rohrlich, would never suspect such refinements are necessary. 

Another attempt to cope with the overlapping divergences was made by Salam 
\195Vo\ 1 195 la\) . I will not describe it, if for no other reason than that I never have 
succeeded in understanding it. Salam and Matthews commenting on this and related 
work somewhat later (1951) remarked "... The difficulty, as in all this work is to find a 
notation which is both concise and intelligible to at least two people of whom one may 
be the author". The belief is widespread that when Salam's work is combined with later 
signihcant work by S. Weinberg il I960)) , the result should be a mathematically consis- 
tent version of renormalization theory. At least that is what one reads in the text book 
of Bjorken and Drell for quantum electrodynamics (Bjor ken and Drell, 1965^ , and in 
the work of R. Johnson (1970) and the lectures of K. Symanzik for meson theories 
(Symanzi k, 1961\ ). So apparently the Matthews-Salam criterion has been satisfied. I 
only wish they had spelled it out a little for the peasants. 

Another foundation of renormalization theory with a rather different starting point 
was put forward by Stueckelberg and Green (1951]). It was refounded and brought to 
a certain stage of completion in the standard text book of Bogoliubov and Shirkov 
(1959 1 ). The mathematical nut that had to be cracked is in the paper of Bogoliubov 
and Parasiuk(1957), (amazingly, not quoted in the English translation of Bogoliubov 
and Shirkov)]j. This paper introduces a systematic combinatorial and analytic scheme 
for overcoming the overlapping divergence problem. This paper is very important for 
later developments. Unfortunately it was found by K. Hepp (1966) that Theorem 4 of 
the paper is false, and that consequently the proof of the main result is incomplete as it 
stands. However Hepp found that Theorem 4 is not essential to derive the main result 
and he could Gil all the gaps. Thus it is appropriate to introduce the initials BPH to 
stand for the renormalization method described in {Bogoliubov and Parasiuk~^1 957) 
and {Hepp, 1966\ ). So far as I know it was the first version of renormalization theory 
on a mathematically sound basis." 



2 see also l |Parasiuk, 1 9601, 
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The rest of the article is also highly recommendable. Reading his article one can 
understand why many workers in the field who had gone through these historical de- 
velopments were wary about renormalization theory [§. A real breakthrough in the 
ease of understanding was supplied by Zimmermann (|1968[ H973a|, 1 1 9 73 . (in partic- 
ular also concerning renormalization of composite operators). I strongly recommend 
Lowenstein's article (|1972j) for a clear exposition of Zimmermann's methods. 

It is the purpose of this chapter to give an overview of the important steps and 
concepts of perturbative renormalization, stating the main results without proofs; for 
the latter the interested reader must consult further literature. Moreover I will only 
discuss the standard approach using expansions in Feynman diagrams, but I would 
like to mention a powerful alternative approach using flow equations based on the 
renormalization group developed by Wilson (|1971«l 1197161 11974p and improved by 
Polchinski (|1984p . a framework which is well suited for proving structural results, see 



e.g. the works of Keller, Kopper and Salmhofcr (1992), (Salmhofer, 1999 1. 

1.3 Perturbative classification of theories 

Concerning the nature of a Feynman diagram the first important concept is its su- 
perficial degree of divergence. Consider a diagram with E external lines, V internal 
vertices, P internal lines and £ independent loops. There are some relations between 
these the first of which is the topological relation 

V + £ — P = number of connected components . (1-3) 

Other relations depend on the precise nature of the vertices in the theory e.g. with 
only a 4-point vertex: 

W = 2P + E. (1.4) 

Then for a connected diagram T with £ loops the superficial degree of divergence 6(T) 
estimates the behavior of the associated integral when all internal momenta are large. 
For (f)^ theory (in d space-time dimensions): 

8{T)=d£-2P (1.5) 

= d--(d-2)E + (d-4)V. (1.6) 

For d — 4 this simplifies to 5(T) = 4 — E , which is negative for E > 4 indicating 
possible overall convergence of the integral in this case. 

Exercise 1.1 Show S(V) = 6 - 2E for 4>% theory. 



This theory is often discussed in the literature because of the comparative simplicity 
of the diagrams involved. In this case 5(F) is negative for E > 3. 

including my first lecturer (in 1966) on QFT, who began his lectures by saying "There are only 
about 3 people in the world who really understand renormalization theory, and I am not one of them!" 

4 A (sub)diagram is connected if all pairs of its vertices are connected by a path following internal 
lines. 
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Fig. 1.3 2-loop diagram in (f>% theory with a nested divergence 

It is clear that considerations of 6(T) alone are not sufficient to establish con- 
vergence of the integral because T could have divergent sub-diagrams as in Fig. 11.31 
However it does lead to the following classification of perturbative QFT: 

• super-renormalizable: theories where only a finite number of diagrams have 5 > 0. 

• renormalizable: there exists Eq < oo such that all diagrams with E > Eq have 
S < 0. 

• non— renormalizable: all E-j>omt diagrams have 6 > for sufficiently high number 
of loops I. 

Remarks: i). Non-renormalizable theories are not necessarily unphysical or devoid 
of predictive power (an interesting example concerns gravity ( |Donoghue, 1995] ) ). They 
can be good effective theories, such as chiral Langrangians discussed in Goltermann's 
lectures. 

ii) Non-perturbative formulations of perturbatively renormalizable theories may 
be trivial in the limit that the ultra-violet cutoff is sent to infinity 0. Examples of such 
theories are thought to be 4> 4 and QED in 4-dimensions, albeit there is at present no 
rigorous proof of this conventional wisdom. 

iii) There are some rigorous non-perturbative constructions of super-renormalizable 
theories e.g. §\ (Glimm and Jaffe, 1968[), Yukawa theory in 2-dimensions ( |Seiler, 1975[ ), 



b\ ( |Feldman, 1975 ). The Schwinger functions in the topologically trivial sector of 
SU(2) Yang-Mills theory in 4-dimensions in a sufficiently small volume have also been 
rigorously constructed (Magncn et al. 1992, 1993). For all of these cases the structural 
information on renormalization found in perturbation theory carry over to the non- 
perturbative framework. 

In the following we restrict attention to renormalizable theories. A generic pro- 
cedure to put perturbative renormalization on a firm mathematical basis consists of 
three steps: 

1. Introduce an UV regularization. 

2. Construct a mapping of a Feynman integral to an absolutely convergent integral 
when the cutoff is removed. 



5 the usual perturbative assumption that a finite coupling exists in the continuum limit is not 
fulfilled 
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Fig. 1.4 2-loop diagram with overlapping divergences 

3. Show that the map in step 2 is equivalent to a renormalization of the bare 
parameters and fields in the original Lagrangian, which formally ensures the desired 
axiomatic properties of the resulting amplitudes. 

There are many perturbative regularization procedures 0, each with their own 
advantages and disadvantages. In all cases known they give equivalent physical results. 



1.4 Zimmermann's forest formula 

Let us first restrict attention to theories with massive bare propagators. Returning to 
the simple 1-loop example the integral can obviously be made convergent if the 

integrand is subtracted at external momentum p = 0: 



R = 9lJ 



A 

d 4 fc 



1 1 



(k 2 + m 2 )([(k + p) 2 + m 2 ) (k 2 + m 2 ) 2 



(1.7) 



We note that since the subtraction is a constant (independent of the external momenta) 
it can be absorbed in a renormalization of the bare coupling go- 

The problem now is to find the generalization of the subtraction in the simple 
case above for an arbitrary diagram. For this purpose we can restrict attention to one 
particle irreducible (1PI) diagrams, diagrams which remain connected after cutting 
any internal line, since an arbitrary diagram can be constructed from 1PI parts joined 
by propagators. 

Here we consider mappings where subtractions are made directly on the integrands, 
as adopted by Reisz for the lattice regularization 0. Let Jp (p, k) be the original inte- 
grand of a diagram with external momenta p = {pi , P2 , • ■ • } and internal momenta 
k = {ki, ki, . . . }. The regularized integral will then be given by 

Rr(p) = J l[d% R T (p,k), (1.8) 
j 

with 

6 Some involving exotic mathematics such as quantum groups, p-adic numbers,... 
7 In the case of dimensional regularization the subtraction is simpler since the subtractions are in 
this case e— pole parts of sub-integrals. 
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Rr(p, k) = Ir(p, k) — subtractions (1.9) 
= (l-# r >)S r (p,fc), (1.10) 



where is the Taylor operator of degree 8 

''dtf F{p ' )lp -° + ""' «>!'" <;//'," ">/ 



(i.ii) 



Dyson's original proposal f|1949&[) was 

S r (p,fc)= I] (l-4 7) )i"r(p,fe) 

where p 7 are the external momenta of 7. On the rhs of this equation the term involving 
71 is to the right of that for 72 if 71 is nested in 72 i.e. 71 C 72 ■ If sets of lines of 71 and 72 
do not intersect then the order is irrelevant. The proposal cures the problem of nested 
divergences, however no prescription is given how to order overlapping divergences, 
that is divergent sub-diagrams which have non-trivial intersection but are not nested 
such as in Fig. [LllFl . 

The correct prescription presented by Zimmermann (1968) is obtained by expand- 
ing Dyson's product and dropping all terms containing products of terms involving 
overlapping sub-diagrams. The result is called Zimmermann's forest formula: 

Rrip, fc) = E II (-# 7) ) ^(p, fc ) (1-12) 

where the sum is over all forests: 

^ = { 7 | 7 lPI,<f(7)>0; 

7i , 72 € T => 71 C 72 or 72 C 71 or 71 , 72 non — overlapping} . 

The empty set <j> is included as a special forest. 



Exercise 1.2 Give the set of forests for the diagram in Fig. 11.51 (it has 8 elements). 



The following rules are used to specify the momenta flow in T and its subgraphs. 
The momenta flowing in the er'th internal line L a ba joining vertices a, b is specified 

as laba = Paba{p) + K aba (k) with K a ba = J2i £ abaih; £ a bai = 1 if laba € loop Ci\ 

tabai = — 1 if haa € loop Cjj e a bai = otherwise. The momentum P a ba is solved by 
Kirchoff's laws: J2ba p aba + Pa =0. Y,L ab ^ec z r aba Paba = 0. Here "resistances" r a ba 
can be chosen for convenience (including zero) but no closed loop of zero resistances. 

Similarly for subdiagrams 7: V aha — P^baiv 1 ) + ^2ba(^)- Fire! solve for Pj 6(T keep- 
ing same resistances and momentum conservation at the vertex. Then determine the 

K 2ba( k ) fr0m l lba = l aba- 

8 It is amusing to hear Dyson's account of his realization of the problem at http:webstories.com 
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Fig. 1.5 3-loop diagram in <j>% theory 
After applying the forest formula the resulting integrand has a form 

V(k,p, m) 



F(p,k) = 



niLi(^ 2 +™ 2 ) 



where 



h = Ki{k) + Pi(p) . 



(1.13) 



(1.14) 



To proceed it is first necessary to give conditions on such an integrand required to 
guarantee the convergence of the corresponding integral 



F{p) = J d d k 1 ...d d k e F(p,k). 



(1.15) 



Here and in the following we will assume that the internal momenta have been selected 
in a "natural way" so that Ki = fe« for i = 1, . . . ,£ (this is always possible). Consider 
the set 

C = {K u ...,K n }, (1.16) 

and let u\, . . . , u r , V\, . . . , vg~ r r > 1 , be I linearly independent elements of L. Let 
H be the (Zimmermann) subspace spanned by the Ui. Then the upper degree of a 
function f(u,v) with respect to H is given by 



deg u |„/(it,«) = v if lim A v f(Xu, v) ^ 0, oo . 

A— >oo 

and for the integral (|1.15l) the upper degree wrt H is defined by 
deg H F(p) = dr + deg ulv F(p, k(u,v,p)) . 



(1.17) 



(1.18) 



The important statement is Dyson's power counting theorem: F{p) in (|1.15[) is 
absolutely convergent if deg H F(p) < for all subspaces H . 

The theorem was first proven by Weinberg (|1960p . and later a simpler proof was 
given by Hahn and Zimmermann (1968). 

It can be shown that the forest formula yields integrals F(p) — Rr{p) whose 
integrands satisfy the conditions of the theorem. Once this is done it remains to prove 
step 3 of the renormalization procedure (see sect. 1.3). This is not at all obvious 
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from the forest formula as it stands. For this purpose it is advantageous to show that 
the forest formula is a solution to the Bogoliubov-Parasiuk-Hepp recursion relations 
( Bogoliubov and Shirkov 1959 Bogoliubov and Parasiuk 1957 |Hepp 1966] ) 3 



Sr = EIlH 5(7) ) S 

if) 7£-0 



7 ! 



(1.19) 



where ip are sets of disjoint 1PI subgraphs. This formula serves as a basis for an 
inductive proof of step 3; again the proof is too long to present here 0; a tactic to 
prove multiplicative renormalizability of gauge theories will be given in sect. 1.6. 



1.5 Renormalization of theories with massless propagators 

For theories with massless bare propagators one must control possible infra-red (IR) 
divergences. To this end for a Zimmermann subspace one defines a lower degree: 

deg ,J(u,v)= E if limA-^/(Au,v)^0,oc, (1.20) 

u \ u A— fO 

and for the integral (jl.151) : 

deg H ,F(p) = d(£ ~ r) + deg vlu F(p,k(u,v,p)) . (1.21) 

The appropriately modified power counting theorem then states that F(p) is absolutely 
convergent if deg H F(p) < and deg HI F(p) > for all H and for all non-exceptional 
external momenta p (that is a set for which no partial sum of external momenta vanish, 
except for the complete sum expressing total momentum conservation). 

To achieve this one has to modify the i?-operation described above to avoid IR di- 
vergences introduced by Taylor subtractions at p = 0. One way, introduced by Lowen- 
stein and Zimmermann (| 1975[) and adopted by Reisz (|1988al [1988c) is to introduce an 
auxiliary mass in the massless propagators l 2 —> I 2 + (1 — s)p 2 . UV subtractions are 
then made at p = 0, s = 0, and afterwards s is set to 1 in the remaining parts. In order 
to avoid further IR singularities it is usually necessary to make extra finite subtractions 
to ensure 2- and 3-point functions vanish at exceptional external momenta singular- 
ities. In gauge theories this is guaranteed by the Slavnov-Taylor identities which we 
consider in the next section. 



1.6 Renormalization of 4-d Yang-Mills theory 

In this section we will outline the proof of multiplicative renormalizability of Yang- 
Mills (YM) theory in 4 dimensions. Perturbation theory is an expansion around a 
saddle point. For gauge theories, such as pure YM theory, gauge invariance of the 
action implies severe degeneracy of the saddle point. To apply perturbation theory 
one needs to lift the degeneracy by gauge fixing (see Liischer's elegant discussion 

9 It has been noted that the BPH R— operation has the structure of a Hopf algebra; see e.g. Kreimer's 
book (1200011 

10 For slightly simplified proofs of BPH renormalization sec (Caswell and Kennedy, 1982). 
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It is convenient to employ linear gauge fixing functions, e.g. the Lorentz 
gauge, which transform under the adjoint representation of the global gauge group so 
that the gauge fixed action is invariant under such transformations. Because of the 
non-linearity of the YM theory this way of fixing the gauge introduces extra terms 
in the functional integral, which can be expressed as a local Lagrangian involving 
Faddeev-Popov ghost fields c, c as introduced in Hernandez' lectures. The action in 
this case is 

So = -2 |tr|iF 2 + ^(9A) 2 -c^[^,c]| (1.22) 

= S A + S gt + S FP , (1.23) 

where = + ,g ^4 M and F^ u = d^A v -d„A fl + g [A p , A v ] . Here A^ = T d A^ where 
Td denote a basis of the Lie algebra of SU(iV) with tr(T<jT e ) = — hSde- 

The action is no longer gauge invariant but an important property is that Sa and 
Sgf + Spp are separately invariant under BRST transformations (Becc hi et al., 1976[ ): 

8A tl = esA„ = e[D tl ,c], (1.24) 

5c = esc = —egoc 2 , (1.25) 

5c = esc = eXod^Aft , (1-26) 
where e is an infinitesimal Grassmann parameter. 



Exercise 1.3 Show the BRST irivaxicince of .S'gf -t- iSfp • 



Note that s is nilpotent: s 2 H[A,c,c] — for all functionals H. This property and 
the BRST invariance of the action and measure imply the powerful Slavnov-Taylor 
identities ( |Slavnov, 1972[ ), flTaylor, 1971] P| for any functional H[A,c,c\. 



For example if we take H 



(sH) = . 
c a (x)d v A h v (y) we obtain 



{\ d^Al{x)d v Al{y)) = {c a {x)d v [D u ,c] b (y)) 
= 5 ab 5(x - y) , 



(1.27) 

(1.28) 
(1.29) 



where in the last line the ghost field equation of motion has been used. It follows that 
the full (bare) gluon propagator has the form 



5ab 
k 2 



i + n(F) ' A k 2 



1 h^h v 



(1.30) 



More generally BRST invariance gives relations between bare Green functions and 
those containing composite operator insertions which describe the non-linear symmetry 

11 For gauge invariant correlation functions it is in principle possible to do PT in finite volumes 
without gauge fixing by solving systems of coupled Schwinger-Dyson equations. But in practice this 
is not usually done because the system is difficult to solve analytically (without hints of the structure 
of the associated Feynman diagrams) . 

12 According to Muta H1987I I. these generalized WI were first discussed by 't Hooft (1971). In his 
discussion only a part of the system of identities was taken into account 
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transformations of the fields. Many features of the theory derive from these Ward 
identities. In particular they are sufficient to restrict the structure of the counter- 
terms and show multiplicative renormalizability i.e. the theory is renormalized by 
adjusting the bare parameters and fields of the bare action. Also the renormalized 
theory satisfies similar symmetry properties as the bare one lH, which is also crucial 
to prove the IR finiteness of correlation functions with non-exceptional kinematics. 

To proceed systematically the WI's are summarized in equations for generating 
functionals ( Becchi et al., 1976] ), ( |Zinn-Justin, 2002] ) such as 

Zo{J,Z,t,K,L;g ,\o)=N J D[A]D[c]D[c}exp{£ (A,c,c,K,L;g , A ) 

+ ~S c (A,c,c; J, (1-31) 

where No is chosen such that Zo(&\\ sources = 0;go,Xo) = 1- The term S c contains 
source terms for the basic fields 

S c (A,c,i»J,lt) = f ^,{4(x)Al(x)+e (x)c b (x)+c b (x^ b (x)} , (1.32) 

Jx b 

and 

8 (A,c,c,K,L;g ,X ) = ^ f £ {K^sA^x) + L b (x)sc b (x)} 

~t s o(A, c, c;g , A ) , (1.33) 
n 

where the new sources K, L must be introduced for the composite fields appearing in 
the BRS transformations. Note that since s is nilpotent further sources in Zq are not 
needed. The factors of fi in diagrammatic contributions to correlation functions count 
the number of associated loops. 

Now the BRS transformations imply that Zq satisfies 



&lJ^-*oe(x)d fll r£ rT }Z =0, (1.34) 



dK»Jx) ^ b dL\x) v '^dJ b (x) 



and the same equation holds for the generating functional of the connected Green 
functions Wq = TivtiZq. For considerations of renormalizability we have already noted 
that it is more convenient to consider 1PI diagrams. The generating functional Tq for 
these is called the vertex functional (see Hernandez' lectures) which is obtained from 
Wq via a Legendre transformation, 

Wo(J,£,€,K,L',g ,\o) = T (A,c,c,K,L;g ,Xo) 
+ I Wx)Al(x) + e(x)c b (x) + c b (x)e(x)} , (1.35) 

J X 

13 since the symmetry transformations are multiplicatively renormalized, 
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where 

A b (x)- dW ° c b (x)- 9W ° c b (x) - dWo (136) 
with inverse relations 

J » = , ?(*) = 7^ > = -J^Ts ■ (1-37) 

The vertex functional satisfies the equation 

f dT dTp dT dT b dT \ 

Setting 

r = r + / (Ao/2) (d^Or)) 2 , (1.39) 

and using the FP ghost field equation 

dT dT , . 

^M^ + c^T°' (L40) 

we can write the WI in a slightly simpler form 



( dV dTp dV dV Q 
\dA b (x) dK b {x) dc b (x) dL b (x) 



= 0. (1.41) 



We now seek a functional £r = £n + 0(h) such that Green functions generated 
from 

Z n (J, C, i; .90, A ) = iV R /" D[A]D[c]D[c] exp{£ R (A, c, c, if, L; g, A) 

+(l//i)5 c (A,c,c;J,f,0} (1-42) 

are finite order by order in PT Such a functional can be found and is simply 
characterized by the following 

Renormalization Theorem: There exist renormalization constants Z\, Z3, Z3 which are 
formal power series in h, relating bare fields and renormalized fields 

A° = Zl /2 A , L° = Zl /2 L , A = Z- x \ , g Q = Z l9 , (1.43) 
c° = Zl /2 c , c° = Zl ,2 c , K° = Zl /2 K , (1.44) 

such that 

£ R (A, c, c, if, £; 50, A ) = £ (A°, c°, c°, if , i°; 3o , A ) . (1.45) 
Furthermore the Z's have only (multi-) poles in e with dimensional regularization (DR) 



't Hooft and Veltman, 1972] ), (jBreitenlohner and Maison 1977ai[l977&|) . ( |Collins, 19841 ), 



or depend only logarithmically on a momentum cutoff with other regularizations. 

14 Clearly £r isn't determined uniquely by this condition 
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The renormalization theorem implies Tr(A, c, c, K, L; g, A) = To (A , c°, c°, if , L°; <?o, Ao) 
and the renormalized vertex functional satisfies the same equation as IV The form 
of the WI indicates universality of the renormalized theory, (the correlation functions 
being fixed by a finite number of specified normalization conditions). The measure 
and renormalized action Sb,(A, c, c, g, A) = So(A°, c , c°, go, Ao) are invariant under 
5 e A p = Z^^SeA ^ , . . . The same holds true for the contributions to £r which are 
linear in K, L. 

The proof of the theorem is by induction on the number of loops. Assume the 
theorem holds to order n — 1 in the loop expansion and hence the functional 
satisfies the identities. The key point is that the equations for T are non-linear, and 
one obtains simple equations (involving functional derivatives of the action) for the 
divergences in the limit of large cutoff which occur at order n. Since all subdivergences 
have been removed, only overall divergences remain which by the power counting 
theorem are polynomials in momentum space of the order of the divergence degrees of 
the corresponding diagrams 0- The equations from the WI suffice to restrict all the 
coefficients (i.e. the renormalization constants) such that the theorem can be proven 
at order n. 

From the renormalized BRS WI for the functionals follow various relations for 
correlation functions. Starting from the identities (jl.38ll.41P for To or Tr and differ- 
entiations wrt the FP fields and the gauge fields. As mentioned before one gets e.g. 
Slavnov-Taylor identities for 2-,3-point vertices, which imply absence of mass terms in 
the gluon self-energy and small momenta behavior of 3-point vertex function. 

1.7 Perturbative renormalization of lattice gauge theory 

For an introduction to perturbation theory with the lattice regularization please see 
the lectures of Hernandez. Here I just summarize some of the salient points. First the 
lattice action is chosen to be gauge invariant and such that it has the desired classical 
continuum limit. The Feynman rules are algebraically more complicated than in the 
continuum e.g. one also encounters vertices involving more than 4 gauge fields, but 
they are straightforwardly derived and now usually computer generated. There are 
also extra terms coming from the measure (the Jacobian of the change of variables U 
to A). The integral over the gauge field A is extended to oo, which formally modifies 
the integral by only comparably negligible non-perturbative terms. For performing 
the perturbative expansion in this way, gauge fixing is needed as in the continuum. 
Finally there is an exact BRS symmetry on the lattice since this symmetry doesn't 
originate from the continuum aspects but is a general consequence of gauge fixing in 
the presence of a gauge invariant cutoff (see Liischer (1988)). 

So we have practically all the ingredients to prove perturbative renormalizability 
following the continuum procedure apart from a power counting theorem for the lattice 

15 Degrees of divergence are denned as before. The superficial degree of divergence of a diagram with 
JVq, Npp , Np external gluon, ghost, and fermion lines is: d = 4 — Nq — ^(Npp+Np) have 8 cases with 
d > 0: (N G ,N F p,Np) = (0, 0, 0), (2, 0, 0), (0, 2, 0), (0, 0, 2), (3, 0, 0), (1, 2, 0), (1, 0, 2), (4, 0, 0). The vac- 
uum diagram can be dropped because it is absorbed in the normalization of the generating functional. 
The rest have degrees 2, 1, 1, 1, 1, 0, 0, respectively. Although superficial degrees > the full ampli- 
tudes are all only log divergent due to the gauge symmetry and Lorentz invariance. 
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regularization. The latter has been provided by Reisz (11988 &[) : very nice accounts can 
be found in ( |Luscher, 1988D and QReisz, 1988d| - An influential earlier paper consider- 
ing renormalization of lattice gauge theories is by Sharatchandra (|1978p . The difficulty 
comes from the fact that in momentum space the domain of integration is compact, 
internal momenta k are restricted to the Brillouin zone B, [k^l < 7r/a, and the inte- 
grand is a periodic function wrt the loop momenta. Also, as mentioned above, often 
the vertices are quite complex and there are additional "irrelevant" vertices i.e. those 
which vanish in the naive continuum limit. In this situation the usual convergence 
theorems in the continuum do not apply. 

Let us first consider the case with massive propagators. A general lattice Feynman 
diagram corresponds to an integral of the form Qf| 

I = / d fci . . . d k t — , (1.46) 

Jb C(k,p,m,a) 

where the numerator V of the integrand contains all the vertices and numerators 
of propagators, and the denominator C is the product of the denominators of the 
propagators, 

In order to make rigorous statements on the convergence Reisz specified the fol- 
lowing restrictions on V, C which however hold for a large class of lattice actions: 
For V: 1) 3 u> st V(k,p,m, a) = a~ u F(ak,ap, am) , 

where F is 2n periodic in aki and a polynomial in am. 
2) lim a _>o y(k,p, m, a) — P{k 1 p,m) exists. 

ForC*: 1) C = 11™= l Ci(h,m, a) , with 

Ci(h, m, a) = a~ 2 Gi(ak, am) , (1-47) 
k(k,p) = ajjkj + y] bupi , aijETL, (1-48) 

3 I 

the latter (integer) condition maintaining the 27r-periodicity of G^. 

2) \ira a ^ Ci(l h m,a) = l\ + m 2 - 

3) 3 a , A st \Ci(k,m,a)\ > A{lf + m 2 ) Va<a ,Vi, 

the latter condition ensuring that l/|Cj| ~ 0(a 2 ) at the boundary 
of B. 

With these specifications, for a > and m 2 > the Feynman integrals are abso- 
lutely convergent and the dependence on external momenta is smooth. However the 
lattice integral does not necessarily converge in the continuum limit if the continuum 
limit of the integrand is absolutely integrable e.g. 

I^^I^^ + OCA (1.49) 
(27r) 4 k 2 + m &a z 

i.e. the integral is quadratically divergent although the continuum limit of the inte- 
grand vanishes. 

16 where for simplicity of notation we consider the case of just one mass parameter m. 
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A lattice degree of divergence involves the behavior of the integrand for large 
internal momenta fcj and simultaneously for small a. For a Zimmermann subspace 
(defined as before) the (upper) degree of divergence is defined by 

deg^F = v if F(Xu, v, m, a/A) ~ a^oo k\ u , (1.50) 
degul^ 7 = 4r + deg„|„V-degV| y C'. (1.51) 

With this definition of the degree of divergence Reisz proved the Theorem: If 
&eg H I < for all Zimmermann subspaces, then the continuum limit of / exists and is 
given by the corresponding continuum expression. 

To prove the renormalizability one must also specify a lattice i?-operation corre- 
sponding to subtractions of iocaJ lattice counter-terms e.g. 

t s -> i s , i s p f( P ) = /(o) +p M A/(p')| p , =0 + ... (i.52) 

where fc M = a -1 sin(a/c M ) is 2ir periodic in afc. 

The proof of the theorem is lengthy and cannot be presented here. Let me just men- 
tion that the subtractions are organized by the Zimmermann forest formula as in the 
continuum, and zero mass propagators are treated as by Lowenstein and Zimmermann 
mentioned in sect. 1.5. Some fine points involve the treatment of the measure terms. 
The use of the BRS Ward identities also proceeds similarly to that for continuum 
regularizations QReisz, 1989[ ). 

Finally matter fields are included straightforwardly. Note that criterion 3) is satis- 
fied for Wilson fermions (which require an additional additive mass renormalization) 
for which Ci(p) = (1 + am)p 2 + m 2 + ^a 2 Y^u<uPfiPi>- ^ is not satisfied by staggered 
fermions, however an extension of the power-counting theorem in this case has been 
supplied by Giedt (12007D . This paves the way for a proof of perturbative renormaliz- 
ability of staggered fermions 



17 This remains to be done, but Giedt anticipates no additional principle difficulties (private com- 
munication). 
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Renormalization group equations 



In the last lecture we learned that correlation functions of the basic QCD fields 
are multiplicatively renormalizable provided bare parameters are also multiplicatively 
renormalized. For perturbation theory in the continuum we usually employ dimen- 
sional regularization where Feynman rules are developed in D = 4 — 2e dimensions 
d Collins, 1984[ ). Bare amplitudes have (multi-) poles at e = and in the minimal 
subtraction (MS) renormalization scheme one subtracts just these. Renormalization 
constants then have the form 

Z = l + g 2 ^ +5 4 (| + 7)+--- (2-1) 

where g is the (dimensionless) renormalized coupling related to the bare coupling by 

g 2 = tx- 2e g 2 Z g , (2.2) 

with /x the renormalization scale. Note the MS scheme is a mass independent renor- 
malization scheme i.e. a scheme for which the renormalization constants do not depend 
on the quark masses. 

Renormalized correlation functions involving r gauge fields, n quark- antiquark 
pairs and I ghost -antighost pairs are given by 

G^(M,P;S , ,A,m, J ) = (z^ 1/2 ^j (Z~ 1/2 ) (X~ 1/2 ) G r ' n ' l (p;g , X ,m oj ) (2.3) 

where the renormalized gauge parameter and renormalized quark masses are given by 

A = Z 3 Xq . rrij = Z nl m 0j . (2.4) 

The renormalization group equations ( jCallan, 1970| ), flSymanzik 1970|fl97ip . follow 
immediately from the simple observation that bare correlation functions are indepen- 
dent of the renormalization scale fj,: 

G^ n ' l = 0, (2.5) 
where the coefficient functions are defined through: 



fi— + p— + rrrij- h 5— + 7-73 + 2nj 2 + 2/73 

Ofi og dm* o\ 
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P(e,9) = ^\g ,X ,m 0j = -eg\l- -g—lnZ g 



-eg + /3(g) • 



1 d 



miT(g) = ^ 



di 



d\ 
1 d 

73(5) = o^^T lnZ 3lso,Ao,mo 3 , 



5 (9) = ^—\g M,m 0j , 



and similarly for 72,73- Note that for any function / depending only on g, e, 

d d 
V-Q^f(g,£)\g ,\ ,m 0j = 0(e,g)—f(g,e) . 

The functions /? and r have the following perturbative expansions: 

00 

P{g) = -g 3 J2 b ><9 2k , 

k=0 

r(9) = -5 2 E4 5 2fe , 
with leading coefficients given by 
bo = (4tt)- 2 
bi = (4tt)- 4 
d a = (4tt)- 



fe=0 



3 3 



34 



3 

_ 2 3(iV 2 -l) 



13 



— N z - —N -—)Nf 



N 



N 



(2.6) 
(2.7) 
(2.8) 

(2.9) 

(2.10) 

(2.11) 

(2.12) 
(2.13) 

(2.14) 
(2.15) 
(2.16) 



Note that in the MS scheme f3 and r are independent of the gauge parameter A. To 
show this start from the general relations between renormalized and bare quantities 
in the form: 



g = LL e <^g + ^2a r (g,m k /ii, A)e r | , 

OO 

m j = rrij + fJ,^2 b jr(g, mk/fi, A)e~ r , 

r=l 

00 

Ao = \ + ^2c r (g,m k /n, A)e~ r , 

r=l 

where the coefficients (3i7" . bjf - Ct 1 £11*6 independent of e. Now define 



(2.17) 
(2.18) 
(2.19) 
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_ dg dX _ drrij d\ 



(2.20) 



where the differentiation wrt Ao is at fixed go,moj,p. For any renormalized Green 
function G of gauge invariant operators we must have an equation of the form: 




+ o-g } G = , 



(2.21) 



i.e. a change in the gauge parameter must be compensated by a change of the renor- 
malized parameters and a multiplicative renormalization of the operators appearing 
in the definition of G. Now differentiate eqs. (|2.17|) . (|2.18j) wrt Aq to obtain 



.4 



p 



v = , v = 



_ ( Er=l %X 6 



EOOj 
r=l ~d\~ 



(2.22) 



where 



.4 



1 + J2r=l 9 dc 

Er=l 



dbjr 

dg 



Er=l 



da r 



Ojk 



Edbj 
r=l 



(2.23) 



d(m k /p,) 



Since p and Vj must be finite we conclude that A v must also be finite. It follows 
that since a r , bj r do not depend on e we must have 



p = , v,j ■ = V j . 



(2.24) 



and hence 



da r 
~dX 



dbj r 



OVj 



(2.25) 



The functions /3, Tj are determined by ai,bj% and are hence also independent of A. 

2.1 Physical quantities, A-parameter and RGI masses 

A physical quantity P is independent of wave function renormalization and indepen- 
dent of A and hence satisfies the simplified RG equation: 



^ dp, dg J drrij 



P = Q. 



(2.26) 



Now every solution of this equation can be expressed in terms of special solutions. 
Firstly the A-parameter which doesn't involve quark masses: 



A = pl{g) . 



(2.27) 



where £(g) satisfies the equation 



1 + /3 



dg. 



e(g) - o . 



(2.28) 



Physical quantities, A-parameter and RGI masses 19 



and is completely fixed by its behavior for small g: 



e(g) = ( 6o5 2 ) -6 1 /(2^) e -l/(2 6off 2 ) 



exp 



d.r 



1 1 &i 

P(x) b x 3 b^x 



The other parameters are the RG invariant masses 

Mi = mi 9(g) 

where 9(g) satisfies 

^ + T %) = o, 

which is also fixed (i.e. also its normalization) by its behavior as g — > 0: 

~t(x) d 



6(g) = (2b g 2 )- d °^ 



exp 



d.r 



j3(x) b^x 



(2.29) 

(2.30) 
(2.31) 

(2.32) 



Consider as an example the case when P(Q 2 , u., g,mj) is a dimensionless physical 



quantity depending on a Euclidean momentum Q then 

P = P(Q 2 /U 2 ,g,m j /VQ2~ 



(2.33) 



We can now appreciate the power of the RG equation; it enables us to deduce the 
behavior of P for large Q 2 , since it implies 



P = P !,£(*), m^tVVQ 2 , t = \n(Q 2 /A 2 ). 



(2.34) 



where g(t) is a running coupling and rrij(t) is the running mass. The running coupling 
satisfies the equation 

Q^=P(9), (2.35) 



and is implicitly defined by 



t 



bog 2 



+ ^ \n(bog 2 ) + 2 f dx 



1 1 &1 

P(x) box 3 b\x 



As Q 2 — > co the running coupling tends to zero (logarithmically) 



1 



b t 



bi 



b 2 t 



(2.36) 



(2.37) 



a property known as asymptotic freedom QGross and Wilczek, 1973D , flPolitzer, 1973[ ). 
The running mass is given by 



m 3 (t) = M^ 1 (g(t)) , 



(2.38) 



and satisfies the equation 
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Q^-=fn j r(g). (2.39) 
For Q 2 —> oo it decreases according to 

/o \ <W(26o) (Ah d 1 

, %(( ) . M, (?) (l - + ,„„)) + j* + 0(0} ■ P-40) 

Any two mass independent schemes can be related by finite parameter renormaliza- 
tions; 

9'=9\JxM, X S (<7) = 1 + X<V + --- , (2.41) 

m' = m Xm (g)- (2.42) 

Physical quantities are scheme independent and so (for the case where the renormal- 
ization scales are equal) 

P{^g,m ] ) = P'{^g , ,m l J ), (2.43) 

satisfies (|2.26[) in the first scheme and a corresponding equation in the second scheme 
with new coefficients (3' , r' which are related to /3, r by: 

W) = {/%)§^}l g = s C S ')> ( 2 - 44 ) 
r'(g') = |r( 5 ) + m§- g lnXm(5)} \ g = g ( g >) • (2.45) 



Exercise 2.1 Show that it follows that b'o = &o,&'i = &i,do — <^o i.e. these coefficients are 
universal but the higher loop coefficients are not e.g. 

d[ =di + 2b oX £ ) - d oX f ] • (2.46) 

Also 

A '= Aexp {€}' (2 - 47) 

M' =M , (2.48) 

i.e. A-parameters are scheme dependent (albeit their relations require just 1-loop computa- 
tions), but the RG invariant masses are scheme independent. 



2.2 The MS lattice coupling 

With the lattice regularization of QCD we can also define a MS scheme by just sub- 
tracting powers of logarithms in the lattice cutoff a. The MS perturbatively renormal- 
ized lattice coupling then has the form 

.giatt =go- b Q gl ln(a/i) + . . . (2.49) 

The 2-loop relation between gutt and g^g has been computed for a large class of 
actions, also including fermions using sophisticated algebraic programs for automatic 



Renormalization of composite operators 21 



generation of the Feynman rules (see e.g. QLuscher and Weisz, 1995D , flHart et al, 2009D , 
flConstantinou and Panagopoulos, 2008D ). As discussed in the last section, the 1-loop 
relation gives the ratio of A-parameters. This was first computed for the standard Wil- 
son action by Hasenfratz and Hasenfratz (1980); in that case for N = 3,N{ = one 
finds a large number Aj^g/Ai at t = 28.8. Close to the continuum limit, assuming that 
<?o = is the critical point which should be approached when taking the continuum 
limit, we should find for a physical mass m in the limit of zero mass quarks 

m/A att = c m + 0(aP) , (2.50) 
4ma~^ cM)^ l/a 'e- 1/(24 » 9 ^(!lo), R(g ) = 1 + O(g 2 ) . (2.51) 

An observation of the leading behavior above is called "asymptotic scaling" ; present 
spectral measurements are considered consistent with these expectations although the 
correction factor R(go) is a power series in go and hence usually not slowly varying in 
the regions where the simulations are performed. 

We remark here again that for lattice regularizations which break chiral symmetry, 
such as Wilson fermions (see Hernandez' lectures), the quark mass needs also an 
additive renormalization 

to = Z m m q , m q — too — m c ■ (2.52) 

Computation of the quark self energy to 1-loop gives am c , Z 2} Z m to 0(<7q), and sub- 
sequent comparison to the result in the MS scheme yields Xm'- 

2.3 Renormalization of composite operators 

Consider a perturbative computation of a correlation function involving a bare com- 
posite operator and a product of basic fields at physically separated space-time points. 
After performing the required renormalization of the bare parameters and basic fields 
the resulting expression either 1) stays finite as the UV cutoff is removed or 2) di- 
verges. If case 1) holds to all orders of PT then this is (usually) due to the fact that 
the composite operator is a conserved current or one satisfying a current algebra. In 
case 2) the simplest situation is that the correlation function becomes finite when the 
composite operator is multiplicatively renormalized, 

<h.(jA) = ^(M)0 baIC • 

But in general operators mix with other operators having the same conserved quantum 
numbers and the same canonical dimension (or less) 

Ra = z <yrO° T &xc + "Z x lower dimension ops" . (2.53) 

T 

A correlation function involving an insertion of a purely multiplicatively renormalizable 
operator with a product of diagonally multiplicatively (gauge invariant) renormalized 
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operators, all located at physically separated points, satisfies the RG equation (for the 
case of massless quarks) : 



where 



I4>i 



7r 



// 



(2.54) 



(2.55) 



For the physical interpretation it is often advantageous to define RG invariant opera- 
tors by (in the simpler cases) 



4>-rgh = Ci(/x/A)0 Ri 



where C\ is a solution to the equation 



— +(3— 



Ci=0. 



It is given (with conventional normalization) by 



d(A*/A) = (26o.9 2 (A i ))^ >/(2bQ) 



exp ■ 



where we have assumed that 7^(5) = 7^ 5 2 + 0(.g 4 ). 

As an example of a lattice regularization, let us consider the case of Wilson's 
fcrmions. It has an isovector current 



dx 



70, (x) 



-Y (0) 



P(x) b x 



(2.56) 



(2.57) 



(2.58) 



W = ±{#r)y( 7M - l)tf(s,/i)#r + a/i) 
+#r + afi) y( 7AJ + l)U{x,u)^{x)} 



(2.59) 



where r a are the Pauli matrices acting on the flavor indices, which is exactly conserved 
(in the case of degenerate quark masses) i.e. on shell, d*V"(x) = where 9* is the 
lattice backward derivative. It follows that this bare operator doesn't need any renor- 
malization. No analogous conserved axial vector current exists for Wilson fcrmions 
even for m q = 0. Often in practical numerical computations simpler lattice currents 
are employed: 



V«{x) = i>{x)—^{x) , A$(x) = ^( x )— ltll ^{x) , 



(2.60) 



which are expected to be conserved up to lattice artifacts. In PT this is indeed the case, 
but in order that the currents obey the correct current algebra (up to 0(a) artifacts) 
they require a finite renormalization 



V R = Z V V, A R = Z A A, Z v /Z A = 1 + O{g 2 ). 



(2.61) 
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Digression: It is probably not so well known among students specializing in lattice 
theory that there is a problem with "naive 75" in dimensional regularization. Namely 
the algebraic rules 

{7m'7iJ = 2( W> <W = D > {7m>75} = 0, tr7 5 =0, and cyclicity of trace 

imply the unwanted relation tr (757^7^7,3 7a) = unless D = 2 or D = 4. The 
modified algebra proposed by 't Hooft and Veltman (|1972|) is to define 75 = 70717273 
so that 

{in, 75} = for u, < 3 , [jfj,, 75] = for \i > 3 . 

The algebra is then consistent but results in the necessity of having to introduce an 
infinite renormalization for the bare axial current 

1 1 AT 2 — 1 

^A S (9) = l+9 4 ^26oC F - + ..., C F = ^^, (2.63) 

and we need a further renormalization xa(.<?) = 1 — g 2 ra^ps Cf + - • ■ to obtain a correctly 
normalized current {^l-^(xaZ^ s ) = 0). 

2.4 Ward identities 

A way to renormalize currents, especially also non-perturbatively, is to enforce the 
continuum Ward identities ( Ward, 1950[ ) , ( |Takahashi, 1957| (in some cases only up to 



lattice artifacts), which for the case of flavor SU(iVf)xSU(iVf) are equivalent to current 
algebra in Minkowski space. General Ward identities are obtained by making infinites- 
imal transformations in the functional integral. For transformations which leave the 
functional measure invariant we obtain relations of the form 

(5SO) = (SO). (2.64) 

For axial transformations (e.g. for N{ — 2) ( jLuscher et al, 1996| ) 

6 A ^(x) = oj a (x)^r a l5l p(x) , &j$(x) = w a (x)^(x)^r a l5 , (2.65) 

and working formally in the continuum (with an assumed chiral invariant regulariza- 
tion) the action transforms as 

S A S= [ d i xuj a (x)[-d t2 Al(x)+2mP a (x)] , (2.66) 
Jn 

where P a = ip 1 ^-"/^ is the pseudoscalar density and we have assumed that w a (x) = 
for x outside a bounded region 7Z. For example if the observable O = O cx t has no 
support in 1Z then the WI becomes 

([-d^(x) + 2mP a (x)] Ocxt) =0, (2.67) 

the famous PC AC relation which has many applications (see sect. 5.3). 

lr To prove this start we use the rules to obtain -Dtr757 M 7„ = ^757^7^7^ = — tr7p757 P 7 M 7„ = 
— tr757 / ,7 /i 7„7 P = (— D + 4)^757^7^ so that tr757 p 7„ = unless D = 2. Next repeat the manipula- 
tions with the trace involving 75 and 4 gamma matrices. 
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Fig. 2.1 1-loop diagram contributing to a ipTij) vertex function 
If O — OintOcxti having support inside and outside 7Z, we obtain 

/ d 4 xuj a (x)([-d^(x) + 2mP a (x)] O iDt O cxt ) = (S A O int O ext ) , (2.68) 

which in the limit oj a (x) — > constant inside 7Z, and m = simplifies to 

io a [ da^(x) (A«(x)O int O e rt) = -(S A O int O cxt ) , (2.69) 

J dlZ 

where dcr^ is the outward normal to the boundary of TZ. For the case 0j nt = A h v (y) 
and TZ the region between two fixed time hyperplanes at t 2 , t\, (and using SaA^(x) — 
-icu a (x)e abc V^{x)) the WI reads (t 2 > 

J d 3 x([A a ( x ,t 2 ) - AS(x,ti)]4(y)e>ext) = ie abc (V„ c (y)O cxt } , (2.70) 

which is equivalent to the current algebra relation [Aq(x, t), Aq(y, f)l = i6( 3 '(x — 
y)e abc V c {x,t) in Minkowski space. All CA relations can be obtained analogously. 



2.5 Scale dependent renormalization 

Most composite operators require scale dependent renormalization. The non-perturbative 
renormalization of such operators will be postponed to later sections. Here we just 
outline the simple perturbative 1-loop computation of the anomalous dimensions of 
bilinear quark operators tj)Tip (e.g. the pseudoscalar density T = 75 or V = r°/2), in 
the continuum MS scheme. At 1-loop the bare vertex function of the density with a 
quark-antiquark pair is given by the diagram with a gluon exchange in Fig. 12.11 using 
dimensional regularization 



T(p,p') = -T a T a g 2 J 



d D k 7„7(/s + p)Tj(k - p')7 Jll 1 



(2ir) D (k+p) 2 (k+p') 2 k 2 



A.- 2 

(2.71) 

Noting that the integral is only logarithmically divergent for D = 4, for the compu- 
tation of the UV divergent parts we can set p = p' = in the numerator. Without 
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any previous experience of dimensional regularization one can accept that the integral 
involved is singular as e — > 0, 



d D k C 1 



and it immediately follows for the integrals appearing in (I2.71[) . 

iDu U L (J 



1 



(47T) 



d k k k 



F(k,p,p r ) 



d k k ^kijkpkf- 
(2^rp (fc 2 ) 2 

Then we easily deduce 

r( P ,p')div~-T a 7V 



c 



D(D + 2)e 



(SfivSpT + 2 perms) 



C 



D 



i^ipTipiv - (i - A *)r 



For example d, noting for our normalization T a T a = — Cf 



for T = 7^75 : T(p,p') div = g 



, CC-c 



c 

i Cf 
, CCf 



(2-£) 2 



-(1-A- 1 ) 



A- 1 



7/^75 • 



for T = 75 : rW)div = g 2 ~l5 [d-(i- a- 1 )] 

e 

3 2 77i(3 + A- 1 )75. 



(4tt)V 



(2.72) 

(2.73) 
(2.74) 

(2.75) 



(2.76) 



(2.77) 



To complete the computation of the renormalization factor we have to compute the 
quark field renormalization factor Z2 to 1-loop. I leave this as an exercise; but if we 
accept that the axial current doesn't need any divergent renormalization to one loop 
this contribution must just cancel the contribution above Z^ = ~i and we 

deduce that the divergent part of the pseudoscalar density is given by 



3C F 1 2 



(2.78) 



Note Z, 



(i) 



-do/ 2 as should be the case for ZpZ m = 1. For lattice regularization we 



would obtain Zp = 1 + .95^0 m(a/i) + . . . (in the MS scheme). 



2.6 Anomalies 

Any account on renormalization would not be complete without mentioning anomalies. 
These involve symmetries present in the classical theory which are violated in the 

2 problems with naive 75 algebra do not come in at this stage 
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process of regularization and which cannot be regained in the limit that the UV cutoff 
is removed. This vast and important subject will be covered in the lectures by Kaplan 
and thus no details will be given here. Let me just mention that we have already 
encountered one example which is that massless QCD breaks scale invariance at the 
quantum level. There is a mass parameter in the renormalized theory and the trace of 
the energy momentum tensor is non-zero 



2ff 



(2.79) 



d^A^ = 2mN[^} + J—N[e^ pX F^F pX 



(2.80) 



This anomaly must be reproduced by a given formulation of lattice fermions in order 
that it can be considered an acceptable regularization of QCD. For the many formula- 
tions available the way that this achieved varies quite considerably. For example with 
Wilson fermions the measure is invariant under infinitesimal axial U(l) transforma- 
tions, but the Wilson term in the action breaks the symmetry and produces the correct 
anomaly in the continuum limit as was first shown by Karsten and Smit (|1981[) . In 
the Ginsparg- Wilson formulation it is the measure which is not invariant under the 
(modified) chiral U(l) symmetry transformations (whereas the action is); see sect. 6. 



For a discussion of the energy-momentum tensor in lattice gauge theories see ( |Caracciolo et ai, 1992 \ 

Another famous example (see Kaplan's lectures) is the U(l) axial anomaly ( Adler, 1969) 
expressing the fact that the U(l) axial current is not conserved in the limit of massless 
quarks: 



2.7 Operator product expansions 

Consider correlation functions involving renormalized local gauge invariant operators 
(A(x)B(0)<f)(yi) . . . (f>{y r )) with the j/j physically separated from x, and from each 
other. In the limit x — > singularities appear which are described by local operators 
O n having the same global symmetries as the formally combined operator AB: 



(2.81) 



where the coefficients C^(x) are c- numbers. This so called Wilson's operator product 
expansion has been shown to hold in some generality in the framework of perturbation 
theory by Wilson and Zimmermann (|1972|) . The relation is structural and thought to 
hold also at the non-perturbative level H. If the fields A, B,O n are multiplicatively 
renormalizable then the coefficients obey the RG equation 



d 



d 



cti(*) = o. 



(2.82) 



The engineering dimension of the coefficients are, in AF theories given by those of the 
operators involved. A famous example is the OPE of vector currents 



3 It has been shown non-perturbativcly in some 2d models e.g. in the 2d Ising field theory 
l |Wu et al, 1976"] ) and also in the massless Thirring model (Wilson, 1970). 
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V,(x)Vj(0) ~* =Q C$(x) + CW(x)7V[F 2 ](0) + . . . (2.83) 

The leading coefficient cjS (x) multiplying the identity operator behaves like l/(x 2 ) 3 
(for x — > 0) up to logs, and (in the case of electromagnetic currents) describes the 
leading high energy behavior in e + e~ annihilation. It has been computed to high 
order in PT (Bai kov et al, 2009] ). Because of gauge invariance the next operators 
occurring here (in the massless theory) have dimension 4 (only one term has been 
exhibited above) and so the corresponding coefficients e.g. Cj?J(x) ~ 1/x 2 . As stressed 
by many authors long ago e.g. by David (|1986j) . this does not mean that there are no 
terms in the associated physical amplitudes behaving like l/{x 2 ) 2 . Unfortunately one 
still encounters the contrary statement in the literature; there is reference to "the 
gluon condensate" as a non-perturbativc effect while forgetting that there can be non- 
perturbative effects in the coefficients C$ (x) and the fact that to my knowledge there 
is at present no regularization independent definition of the gluon condensate! These 
considerations do of course not negate the strength of the OPE; PT gives the leading 
short distance behavior of the coefficients and the sub-leading terms in the OPE give 
the leading effects in processes with non-vacuum external states. 

Another useful example of an application of the OPE is in non-leptonic decays 
in the framework of the Standard Model (e.g. K — > 2tt (Gaill ard and Lee, 1974| ), 
( jAltarelli and Maiani, 1974| ) which will be discussed in detail in the lectures of Lel- 
louch). I would just like to emphasize a few points below and in this discussion neglect 
quark masses. The typical Minkowski space amplitude for initial and final hadronic 
states /, F has the form 



T F/ cx J d 4 xD(x,m w )(F\TJ^(x)J^(0)\I), (2.84) 

involving left-handed currents Jf , J%, where D(x, mw) is the scalar function occurring 
in the W-meson propagator. Since the physical W-meson mass mw is much larger than 
typical strong interaction scales involved, short distances dominate the integral. The 
simplest case to consider is that where the currents involve different flavored quarks; 
in that case the OPE implies 

d 4 xD(x,m w )TJ^(x)J^(0) ~ h a (v/m w ,g)O? suv (0) , (2.85) 
with composite operators 

^rsuv — Orsuv i Orsvu , (2.86) 

Orsuv = N [J^J^] , J£„ = Vv 7m \( l ~ 7s)^ • (2.87) 



The operators renormalize diagonally if one has a regularization preserving chiral 
symmetry 0. Restricting to that case we write the rhs in terms of the RGI operators 
introduced in (|2.56l) as 

4 Because of the difficulty with treating 75 in the framework of dimensional regularization, during 
the rcnormalization procedure one has to include mixing with "evanescent operators" (ones which 
vanish for D = 4). 
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^FK/AJO^fO), (2.88) 

cr=± 

with 

k°(m w /K) = h a (l, g 2 (m w ))/C a (m w /A) (2.89) 
= (2b Q g 2 (m w ))-^^ 2bo [1 + k al f(m w ) + ...]. (2.90) 

The coefficients can be computed in PT e.g. k a i = h a \ — (7^' — j^bi/bo)/ (2bo), where 
h a \ is determined by 1-loop matching of the full amplitude with that of the OPE, and 
7^ 1 ' > is a two-loop anomalous dimension. There remains the important job of the lattice 
community to determine the non-perturbative amplitudes (F\0 RGla (0)\I) . 

Wilson fermions break chiral symmetry and this has the effect, as first pointed out 
by Martinelli (I1984p . that with this regularization the parity even part of ± mixes 
with other operators e.g. 0\ A oc (VV — AA) ± (u O v). Although the number of such 
operators is restricted by CPS symmetry (charge conjugation, parity, and SU(-/Vf) in 
the limit m q = 0), there are still 3 such operators which makes Wilson fermions a 
rather awkward regularization for computing the desired physical amplitudes in this 
case. 



3 

Non-perturbative renormalization 



So far we have mainly considered renormalization in the framework of perturbation 
theory, which for QCD is only applicable to a class of high energy processes. But QCD 
is a candidate theory for the hadronic interactions at all energies. In particular in most 
numerical QCD computations we are attempting to determine low energy observables. 
In such studies we fix the bare quark masses by fixing a sufficient number of scales e.g. 
m-ir/fir , ttik/ fir ■■ to their physical values, and then af 7V (go) gives the lattice spacing 
a(go) for the pion decay constant f T fixed. This is called a hadronic renormalization 
scheme. 

In order to connect a hadronic scheme to a perturbative scheme one could in prin- 
ciple proceed by computing a non-perturbatively defined running coupling, e.g. using 
a 2-current vacuum correlation function, over a large range of energies. Eventually at 
high energies (after taking the continuum limit) we can compare to perturbation the- 
ory and estimate the scales e.g. f^/Aj^. We can proceed similarly for running masses 
and scale dependent renormalization constants. However, despite the huge increase in 
available computational resources and advances in algorithmic development, to mea- 
sure physical high energy E observables with small lattice artifacts and negligible finite 
volume effects still raises the old practical problem that the lattices required need too 
many points (a«l/£«l//,<L). 

Many procedures have been applied in order to attempt to overcome this difficulty. 
The most naive way is to try to use the perturbative relation between the MS coupling 
and the bare lattice coupling ao = 5o/(47r) mentioned in sect. 2.2: 



where the constant k depends on the lattice action. As a non-perturbative input one 
computes a mass scale e.g. a charmonium mass splitting, to give the lattice spacing 
a(ao) in physical units. Now one can use the relation above to obtain an estimate 
for ctMsifJ. = s/a) (for some chosen factor s). This procedure encounters many ba- 
sic problems; firstly it is difficult to separate the lattice (and finite volume) artifacts 
and estimate the systematic errors, and secondly for many actions (e.g. the standard 
plaquette action) one encounters large perturbative coefficients, in fact we have al- 
ready seen this in the large ratio between the lattice and MS A-parameters. It was 
first observed by Parisi (|1985j) Q that large contributions to this ratio come from tad- 
pole diagrams, and that similar diagrams appear in the computation of the average 

1 see also {Martinclli et al, 1982[ ) 



a M§(^) = a + Otldi(afl) + ... 

di (an.) = Sirbo ln(a/x) + k , 



(3.1) 
(3.2) 
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plaquette 

P=±(trU(p))=l + O(g 2 ). (3.3) 

Mean Geld improved bare PT is an expansion in an alternative bare coupling 

a P = a /P. (3.4) 

If one now re-expresses a^g in terms of ap one usually finds that the perturbative 
coefficients are reduced; e.g. choosing the scale /x such that the 1-loop term is absent; 
for the standard action N — 3,N{ — one obtains 

a M§ = a P + 2.f85ctp + . . . , for fia — 2.6 , (3.5) 

with a reasonably small 2-loop coefficient. Computation of a mass scale and the pla- 
quette expectation value gives an estimate of the MS coupling at \i — 2.6/a. There 
are obviously many variants and the technique has been perfected by Lepage and 
Mackenzie (|T993). But I think it is still true to say that systematic errors in these 
determinations are difficult to estimate and the scales achieved are not so high that 
one can be confident to use the result as initial conditions for running the RG equation 
with perturbative beta-functions to higher energies. 



3.1 Intermediate regularization independent momentum scheme 

Another approach is to use an intermediate regularization independent momentum 
scheme which will be mentioned again in sect. 5.3. In this approach one usually con- 
siders correlation functions involving some basic fields and thus one has to fix a gauge 
(and tackle the Gribov ambiguity problem if a covariant gauge is used). Using sim- 
ilar ideas non-perturbative running couplings can be defined as they are in PT e.g. 
from the 3-point gluon vertex function (Alles et at, 1997). For a covariant gauge the 



full gluon propagator in the continuum has the form (II.30[) . and the 3-gluon vertex 
function has the structure (pi + p 2 + P3 =0) 

K\%7M>Pi>P3) = ^f aia2a3 F f , 1 ^ l3 (pi,P2,P3) + ^d a ^D^ l3 (p 1 ,p 2 ,p 3 ), (3.6) 

where /, d are the SU(A^) invariant tensors. At the "symmetric point" (SP) pf — 
M 2 ,PiPj = — \M 2 ,i j we have for the first amplitude 

^iM2M3 Isp = F(M 2 ) {<5 MlM2 (pi - p 2 ) fi3 + 2 perms} + i? MlAI2At3 (3.7) 

where n^R^^^ = if npi = for all pi. The symmetric point was also introduced 
by Lee and Zinn- Justin (|1972p as a set of symmetric non-exceptional momenta in order 
to avoid potential IR problems (at least in the framework of PT). One can now define 
a renormalized running coupling by (jCelmaster and Gonsalves I979fc[ 11979a)) 

ff MO M(Af 2 ) = F(M 2 ) [1 + II(Af 2 )] " 3/2 , (3.8) 



where LI(fc 2 ) is the dynamical function appearing in the full gluon propagator (j 1 . 30|) . 
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With an analogous construction for the lattice regularization, one then measures 
<7mom at various values of allowed M 2 = k f^r) data given value of the bare coupling 
go at which one also measures a mass scale in order to specify M 2 in physical units. 
One then, if possible, repeats the procedure at other values of go in order to attempt 
a continuum limit extrapolation. Then at the largest values of M 2 one has to resort 
to perturbative evolution in order to reach high energies. The method avoids the use 
of bare PT, but for presently feasible lattice sizes one cannot really reach sufficiently 
large M 2 , and in the most optimistic case one only has a small window of M 2 without 
too severe lattice artifacts. 

3.2 Recursive finite size technique 

Lattice simulations are necessarily performed at finite volumes and these effects arc 
a source of systematic error if the measured physical quantities are desired in infinite 
volume. On the other hand we can also make use of the finite volume as a probe of the 
system, as developed in QFT to a high degree by Liischer (1991&, 1991a). For example 
in Aoki's lectures he described how one can extract infinite volume scattering data 
from measurements of finite volume effects on spectra at large volumes. 

Our interest here is to overcome the renormalization problem mentioned above and 
for this purpose it is useful to define a renormalized coupling depending on the volume 
ofv e -g- one determined in terms of the force between two static quarks 

a 99 -(L)cx {r 2 F q g(r,L)} r=L/2 . (3.9) 

There are infinitely many acceptable choices and at large L their behaviors can be 
completely different. The important feature which characterizes them is that at small 
L (where the spectral properties are vastly different from familiar infinite volume 
spectra) we can use PT to compute them as a power series in ajjg (starting linearly) : 

apv(L) = a-^(p) + a-^sdi) 2 [87r&o m(/xL) + c FV ] + • . • (3.10) 

Let us assume that we have decided on the definition of the finite size coupling. 
The tactic to connect a hadronic scheme to a perturbative one, which goes under the 
name of the recursive finite size technique consists of the following steps. 

1) Set the scale on the lattice with largest physical extent L = L max so that L max is 
known in physical units say ~ 0.5fm, and compute aFv(imax)- 

2) Now perform non-perturbative evolution until aFv(^) is known on a lattice of much 
smaller size say ~ 0.005fm. 

3) Assuming perturbative evolution has apparently set in at the scale reached in step 
2, one continues with perturbative evolution and eventually obtains the A-parameter 
in the FV scheme Apv m physical units. 

4) Relate the coupling to the MS-scheme to 1-loop to obtain the ratio of A-parameters, 

2 Note on a given (L/a) 4 lattice k can only take certain allowed integer values. 

3 The RFST was first submitted by Liischer as an (unaccepted) project proposal for the European 
Monte-Carlo Collaboration (EMC 2 ) formed after an early realization that large scale simulations 
required large collaborations. 



32 Non-perturbative renormalization 




0.1 0.2 „/ T 

Fig. 3.1 Extrapolation of Z/max/ro to the continuum limit 

and hence A^jg in physical units. One can then use perturbative running to compute 
a Ms a * an y HE scale say /i = m\y- 

In the following we will outline some of the steps in more detail. The order that 
some of the steps are carried out in practice is not fixed, but let us start with step 1. 
Let us define 



L max : volume where gpy(L maK ) — fixed value 3.48 , (3-11) 

the precise latter value is of course not essential but chosen after initial test runs to 
correspond to lattice sizes ~ lfm. We would like to determine L max in terms of some 
physical unit. Considering first the pure Yang-Mills theory, a convenient quantity is 
Sommer's scale ( |Sommer, 19 94) r Q defined by rQF q g(r ) = 1.65, which in phenomeno- 
logical heavy quark (charmonium) potential models corresponds to a distance ~ 0.5fm. 
To compute L max /ro we first select a value of the bare coupling go on a large lattice 
(say L/a ~ 48) where one can measure r^/a accurately with negligible finite vol- 
ume effects. At the same go one measures <?| v on smaller lattices L/a = 6, 8, . . .16 
and obtains L max /a by interpolation, and hence (L max /ro)(go). Now the procedure 
is repeated at other values and subsequently the data is extrapolated to the contin- 
uum limit using the theoretically expected form of the artifacts (discussed in the next 
section) as illustrated in Fig. 13.11 

Step 2, measuring the evolution was historically done in the reverse order from the 
description above and the standard notation in the following is suited for this. In the 
continuum limit there is a well defined function a, the step scaling function, relating 
the coupling at one volume to that at double the volume: 
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Fig. 3.2 Extrapolation of E to the continuum limit a la Symanzik 

g* v (2L) = a(g* Y (L)) . (3.12) 

With a lattice regularization this is modified to 

& v (2L) = V($ v (L) i a/L). (3.13) 

One then starts with a convenient small lattice, say L/a = 8, and tunes <?o such that 
3py(i) equals some small value u. At the same go one computes <?p V (2L) = H(u,a/L). 
One then repeats this for a manageable range of a/L as illustrated in Fig. 13. 2\ and 
extrapolates the result to the continuum limit thus obtaining a value for a(u) = u' . 
The whole procedure is then repeated but this time starting with an initial value 
g FV (L) = u' (or a value close to «'). After this has been done many (O(10)) times, one 
ends up with a sequence of points for the continuum step scaling function as illustrated 
in Fig. I5~3l 

u k = a(u k+1 ), u = 3.48, (3.14) 

giving g~py(L) at L — 2~ fc i max . At the small values of L one can check whether a(u) 
is well described by the perturbative expectation 

a(u) = u + (2\n2)b u 2 + ... (3.15) 

and if this is the case one can use the beta function with perturbative coefficients to 
compute Apv m physical units. After that the steps are straightforward; a 1-loop com- 
putation relating apy to a-g^ yields the ratio of A-parameters and hence the desired 
value of the product A^jgro relating the scales of the hadronic (LE) and perturbative 
(HE) renormalization schemes. 



2.100 0.11 

1.695 0.056 

1.430 0.029 

1.243 0.014 
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Fig. 3.3 Results for the continuum $fep scaling function 

It is clear that for the success of the RFS method described above, we need a def- 
inition of the coupling which satisfies the following criteria: a) it should be accurately 
measurable, b) it has preferably small lattice artifacts, and c) it should be relatively 
easily computable in PT. 

3.3 The Schrodinger functional 

After some extensive R&D members of the Alpha Collaboration found that cou- 
plings based on the Schrodinger functional (SF) satisfy the above requirements 
( jLuscher et al., 1992] ). It was further realized that the setup is also well suited for the 
computation of renormalization constants in general, and that it is easily extended to 
include fermions and to compute their running masses. 

In this framework one studies the system in a cylindrical volume A with Dirichlet 
boundary conditions in one (the temporal) direction and periodic be in the other 
(spatial) directions (illustrated in Fig. 13. 4p : 

A fe (x,0) = C fe , A k (x,T) = C' k , A k (x + L,t) = A k (x,t). (3.16) 

Formally in the continuum the SF is given by the functional integral 

Z{C,C) = [ D[A]e- s , (3.17) 

JSF be 

4 It was for tunate that M. Liischer was already informed about the SF in scalar theories 
JLuscher, 1985[ l, and that U. Wolff had already suggested a similar construction for the 2-d non- 
linear O(n) sigma model. 
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Fig. 3.4 Geometry of the Schrodinger functional 



which is properly defined so that it is equal to the transition amplitude (C"|e~ HT P|C) 
(H the Hamiltonian, P the projector on gauge invariant states) in the Hamiltonian 
formulation. 

The renormalization of the SF in scalar field theories was first studied by Symanzik 
(1981a). He found that apart from the usual renormalization of the bare parameters 
and fields in the bulk one just requires some extra terms on the boundaries, specifically 
spatial integrals over local fields of dimension < 3. Liischer's paper Q1985I) gives a 
clear introduction to the subject. There are no such local gauge invariant operators 
in pure Yang-Mills theory and so the (bare) SF should in this case not need any 
renormalization besides the usual coupling constant renormalization. We remark that 
one of the first papers considering the related topic of the structure of Yang-Mills 
theories in the temporal gauge were by Rossi and Testa ©80a ( 19806)). 

For small bare coupling go the functional integral is dominated by fields around 
the absolute minimum of the action described by some background field B. The SF 
then has a perturbative expansion 



If the boundary fields depend on a parameter r\ then one can define a renormalized 
running coupling as 



In Z{C, C) ~ T(B) = 9 v 2 r a (B) +T 1 (B) + ... 



(3.18) 




(3.19) 



Regularizing the gauge theory on the lattice the Scrodinger functional is an integral 
over all configurations of link matrices in SU(iV): 




(3.20) 



with the Haar measure and the U satisfying periodic boundary conditions in the spatial 
directions and Dirichlet be in the time direction: 
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U(x,k)\ Xo=0 =W(x,k), U(x,k)\ X0=T = W'(x,k), k = 1,2,3, (3.21) 

where ^ 

W(x,k) = 7>expja^ dt C k (x + (1 - i)afc) | , (3.22) 

i.e. the SF is considered as a functional of the continuum fields C, C" and the continuum 
limit a — > is taken with C, C" fixed. 

One can work in principle with any acceptable lattice action, the simplest being 
Wilson's plaquette action 

S=4y>(p)»tr(l-tf(p)). (3-23) 
So p 

where the sum is over all plaquettes p and the weight w{p) — 1 except for those lying 
on the boundary which is chosen w(p) = 1/2 to avoid a classical 0(a) effect. 
Note that the derivative entering the definition of the coupling is 

dr fdS\ , nnAS 

d-,=W>- (3 - 24) 

The expectation value appearing on the rhs involves only "plaquettes" localized on 
the boundary. These are accurately measurable hence satisfying criteria (a) above. 

As for the particular choice of the boundary fields C, C to make the perturba- 
tion expansion well defined we need the following stability condition: if V(x, fi) = 
exp[aB f _ l (x)] is a configuration of least action (with be C,C) then any other gauge 
field with the same action is gauge equivalent to V. Secondly we would like to have 
criterion (b). How to make optimal choices satisfying these demands is not at all ob- 
vious. Again after some experimentation the Alpha Collaboration made the choice of 
abelian be e.g. for SU(3): 



C k (x) = — I 4> 2 J , ^</> Q = 0, 4> a indep. of x, k , (3.25) 




and similarly for C involving elements (/>[. The induced background field is abelian 
and given by (T = L) 

B Q =0, B k (x) = C+(C'-C)x /L, C = C(t}). (3.26) 

Stability has been proven ( Luscher et al, 1992[ ) provided the (f>'s satisfy 

<f)l < 4> 2 < (f>3, \4>a - I < 27T, 

(and similarly for 4>' a ), and provided TL/a 2 is large enough, albeit the bound not being 
very restrictive e.g. TL/a 2 > 2ir 2 for N = 3. 

With this setup the Alpha Collaboration produced measurements of a running 
coupling (in the continuum limit as far as it could be controlled) over a large range 
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1 10 100 ^[Gev] 

Fig. 3.5 SF running coupling aQu) = cjg F (L)/4-K , fj, = l/L (N { = 0) 

of energies d, as depicted in Fig. 13.51 At high energies the running is consistent with 
perturbative expectations, giving convincing numerical support to the (yet unproven) 
conventional wisdom that the critical coupling is g c — and that the continuum 
limit of the lattice theory is asymptotically free [J Contrary to widespread opinion the 
latter property is non-trivial (so far lacking rigorous proof) and some authors have 
questioned its validity ( |Patrascioiu and Seiler, 2000] ) . 



3.4 Inclusion of fermions 

The inclusion of fermions in the SF framework was first considered by Sint (j!994l 
1995). In the continuum it has been argued by Luscher (|2006p that "natural bound- 
ary conditions" involve linear conditions for the fields of lowest dimension. For Dirac 
fermions these take the form 

S^lbdy = (3.27) 

where, in order to obtain a non-trivial propagator, the matrix B must not have maximal 
rank. If one demands invariance under space, parity, time reflections (xq — > T — xq) 
and charge conjugation, one is left with the possibility 

P+ip(x) = lj)(x)P- = Q , £ = 0; P-tp(x) = ~${x)P + = , x =T, (3.28) 

5 This particular SF coupling runs similarly to PT down to low energies, but this is not a universal 
property of non-pcrturbativc running couplings. 

6 Some additional evidence for the existence of non-perturbatively asymptotically free theories 
comes from studies of integrable models in 2d (see e.g. (jBalog and Wcisz, 2004}). 
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(or P + f-> P_) where P± = |(1 ± 70). Sint showed that with these homogeneous 
boundary conditions the SF is renormalizable without the necessity of including any 
extra boundary terms. 

With the lattice regularization where continuity considerations are a priori missing, 
boundary conditions are implicit in the specification of the dynamical fields (those to 
be integrated in the functional integral) and the precise form of the action close to the 
boundary. For Wilson fermions the terms in the action coupling close to the boundary 
e.g. near xo ~ have the form oc ip(a)P + ip(0) + ip(0)P-ip(a). It is thus natural in 
the corresponding lattice SF to declare fields ip(x), ip(x) away from the boundary i.e. 
< xq < T as the dynamical variables and expect that the bc's fl3.28[) are recovered 
dynamically in the continuum limit. Often in the SF literature one sees the equations 

P+^(0, x) = p(x) , ^(0, x)P_ = p(x) , 
P.i>(T, x) = p'(x) , Tb(T, x)P + = p'(x) . 

These are however not to be considered as specifying boundary conditions, but de- 
scribe couplings of sources for the undefined field components near the boundary. For 
example, defining 

£( X )=P_-A-, £( X ) = _ * p (3.29) 
5p{x) dp(x) 

we can consider correlation functions of the form 

(O a A a (x))~ J[dUd^d^}O a A a (x)e- s \ p=fl=p , = p, =0 , (3.30) 

where all sources are set to zero after differentiating and 

0^-^C(y)^ 75 C(z). (3.31) 
y-7- 

In this setting the extra boundary counter-terms appearing in Sint's original paper 

1/2 

amount to a renormalization of the sources £(x) — » £(x). 

An important point is that the SF fermion boundary conditions imply a gap in the 
spectrum of the Dirac operator at least for go small enough. This has the consequence 
that simulations at zero quark mass m q = with the Schrodinger functional are not 
problematic. 

Also an extra option is to impose quasi-periodic boundary conditions in the spatial 
direction of the form 

■ip{x + Lk) =e l6k ip{x), ~$(x + Lk) = e~ iek 4>(x) , (3.32) 

which are equivalent to modifying the covariant derivative to 

(VfcVO(z) = - \e ia6k/L U{x, k)ip(x + ok) - Mx)] , (3.33) 
a L J 

and similarly for V£. Such boundary conditions with various choices of the 9k serve 
as extra probes of the system. 

7 of the form oc J xg _ d 3 x ip{x)P-ip(x) + J xg _ T d 3 x ip{x)Pj r ij){x) 
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Probably in the future, computers will be so powerful that physically large enough 
lattices will be measurable with very small lattice spacing a, such that lattice artifacts 
become numerically irrelevant. Even so the question of the nature of lattice artifacts 
is of theoretical interest. However, at present it is important in practice to gain insight 
in the form of the artifacts in order to make reliable extrapolations of numerical data 
to the continuum limit. 

Usually we make extrapolations of e.g. ratios of masses of the form (jl.2p assuming 
leading artifacts are predominantly polynomial in the lattice spacing. Lattice artifacts 
are non-universal e.g. the exponent p and the coefficient C\i in (|1.2[) depend on the 
lattice action. This can be used in various ways e.g. if we simulate different actions 
and the data for glueball masses looks as in Fig. 14. 1[ this would be a support of the 
(expected) universality of the continuum limit and one could make a constrained joint 
fit. 
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Fig. 4.1 An old plot of glueball masses in pure YM measured with different actions 
(Morningstar and Pcardon, 1996); the results indicate universality. Green points (plaquette 
action) and the others from improved actions. 
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More ambitious ways of using the non-universality involve designing actions with 
larger values of the exponent p, which are called Symanzik improved actions (jSymanzi k 1979[ 
(| 198 1 6[) ) or even constructing perfect actions having in principle p = oo (see sect. 6). 

Again most of our knowledge concerning lattice artifacts comes from studies of 
perturbation theory. Some non-perturbative support for the validity of the structure 
found there comes from investigations in the 1/n-expansion of QFT in 2 dimensions 
(see (jKncchtli et al., 2005[ ), ( |Wolff et al., 2006] ) and references therein), and also from 
many numerical simulations. 



4.1 Free fields 

Let us first consider a free scalar field theory on an infinite 4-dimensional hyper-cubic 
lattice with standard action: 



5 ° = fl4 E \ [d^{x)d^{x) + m 2 ^) 2 ] , (4.1) 

X,fJ, 

where d fl f(x) = [f(x + aft) — f(x)\ /a. The 2-point function is 



Ar + m z 

where = — sin Noting that for small a, 

fc 2 = fc 2 -^ 2 E fc M + °(« 4 )' ( 4 - 3 ) 

we can write [<f>o the corresponding field in the continuum theory) 

G(k) = G cont (k) - a 2 (SfMk)M0)) + ■■■ (4.4) 

with 

sf = -f d 4 x YA d2 M x)d > o{x) ■ (4 ' 5) 

On-shell information is obtained from the (lattice) two point function G(x) when x 
is separated from by a physical distance. Performing the integral over kg we obtain 
the representation 

7^3 e ikx ° ; , (4.6) 
-Tr/a (?iry R(k,m) 

where the energy spectrum e(k, a, m) is given by 

cosh(ae(k, a, m)) - 1 = ^a 2 (k 2 + m 2 ) . (4.7) 

Defining the pole mass by m p = e(0,a,m) we obtain in the continuum limit a — > 0, 
m p fixed: 
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e(k, a, to) 2 = m 2 + k 2 - — T(k, m p ) + 0(a 4 ) ; T(k, m p ) = ^ fc j + k ' (k 2 + 2 ™ 



(4.8) 

The cutoff effects are (for to p = 0) of order > 10% for k > 27r/(5a). One can improve 
this situation by adding 0(a 2 ) terms to the action. This can be done in many ways, 
the simplest possibility being 

S = So + cSi , (4.9) 



the latter involving interaction of next-to-nearest neighbors. The energy spectrum is 
now given by the solution to 

cosh(ae(k, a, m)) — 1 — 2c [cosh(ae(k, a, to)) — l] 2 = -a 2 [ k 2 + to 2 + ca 4 fc 4 

2 v 

(4.11) 

Now for small lattice spacing a the energy spectrum has the form 

e(k, a, m) 2 = to 2 + k 2 + (c - -^j a 2 T(k, m p ) + 0(a 4 ) , (4.12) 

from which we see that the energy is 0(a 2 ) improved if we chose c = j^. 

Note that for the improved action another energy level is present but its real part 
always remains close to the cutoff, and hence it is irrelevant for the continuum limit 

Exercise 4.1 What is Sl B for the action 

a 2 

So + cSi + d ^2 -^d1d )J ,(l>(x)dtd v 4>{x) . 
Show that the energy is improved for c = for arbitrary d. 



4.2 Symanzik's effective action 

Based on low order perturbative computations in various field theories one arrives at 
the following conjecture: In a large class of interacting lattice theories (in particular 
asymptotically free theories) there exists an (Symanzik's) effective continuum action 

Sf = J d d yC x (y), (4.13) 

such that a Green function of products of a multiplicatively renormalizable lattice field 
ip at widely separated points Xi takes the form 

1 A spectral representation exists but energy levels can be complex. 



42 Lattice artifacts 



Z r J 2 {ip(xx) . . . ¥>(£r))latt = (foixi) ■ ■ • ^o(2v))cont 

-a p J d d y (Citpo(xi) . . . ipo(x r )) cont 

r 

+a p ^((fi (x 1 ) . ..(fi(x k )) ■ ■ ■ <po(x r )) cont + ... (4-14) 

k=l 

where (po,(fi are renormalized continuum fields, in particular (pi is a sum of local 
operators of dimension d v + p depending on the specific operator ip and having the 
same lattice quantum numbers as ip. The effective Lagrangian is a sum 

C l =Y J d{9{ti,w)Om{v), (4-15) 

i 

of local operators Oyh of dimension d + p having the symmetries of the lattice ac- 
tion, /x is the renormalization scale e.g. of the dimensional regularization used in the 
continuum. 

Note a) in the integral over y one in general encounters singularities at points 
y = Xk- A subtraction prescription must thus be applied, but the arbitrariness in 
this procedure amounts to a redefinition of <p\. b) The coefficients Cj are, as indi- 
cated, functions of the lattice spacing a, but the dependence is thought to be weak 
(logarithmic) . 

If the conjecture is true then one generically expects 0(a 2 ) artifacts in pure Yang- 
Mills theory and 0(a) effects with Wilson fermions. All present numerical data seems 
consistent with these expectations but until now only a small range of a is available. 

4.3 Logarithmic corrections to 0(a 2 ) lattice artifacts 

In 2d lattice models e.g. the non-linear 0(n) sigma model, which is perturbatively 
asymptotically free, one can simulate lattices with very large correlation lengths (> 
200a). In these theories the expectation is also 0(a 2 ) artifacts. Hence it came as a 
surprise, as mentioned by Hasenfratz in his LATT2001 plenary talk (Hascnfratz, 2002), 
that data on a step scaling function in this model seemed to show O(a) effects as 
illustrated in Fig. 14.21 This was rather unsettling and motivated Balog, Niedermayer 
and myself (|2010[ I2009|) to investigate the logarithmic corrections to the 0(a 2 ) in 
the framework of renormalized perturbation theory. We found that generic artifacts 
in the 0(n) sigma model are of the form a 2 ln s (a 2 ) with s = n/(n — 2). For n = 1 
the exponent is s — 3, and such strong logarithmic corrections to the 0(a 2 ) effects 
can explain the peculiar behavior, and yield good fits of the data for various actions 
QBalog et at, 2010j ). For ri = oo the exponent is s = 1 which is consistent with what 
is found in leading orders of the 1/n expansion ( [Wolff et al, 2006) . 
The steps involved in obtaining the result above are as follows. 

1) Classify operators of dimension 4 (recall for this case d = 2 and p — 2) which 
appear in the Symanzik effective Lagrangian (|4.15l) . 

2) Compute the Ci at tree level (the coefficients normalized such that Ci = cf^ + 
0(g 2 )). Although finally interested only in on-shell observables, it is sometimes con- 
venient to work off shell and compute a sufficient number of correlation functions G' 1 *' 
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Fig. 4.2 Monte Carlo measurements of the 0(3) a-model step scaling function at 
uo = 1.0595 for two lattice actions. The fits shown contain a and a In a terms. Fits of the 
form kia 2 + kic? In a + k^a 4 have unacceptably high x 2 /dof. 

with a product of r basic fields, in the continuum and on the lattice. For the lattice 
Green function, 



where G- rj are continuum correlation functions involving additional insertion of a 
composite field 0^. 

3) The ratios = G- r ' /G^ which characterize the lattice artifacts (but are 
themselves independent of the lattice regularization) obey an RG equation of the form 



where is obtained from the mixing of the Orj to 1-loop (see sect. 2.3). If we have a 
basis where the renormalization is diagonal to one loop — — 26 A i (5ij<j[ 2 + . . . then 
the operator associated to the largest value of A 4 generically dominates the artifacts 
if the corresponding tree level coefficient Cj(0) ^ 0. 

The program should be carried out for lattice actions used for large scale simula- 
tions of QCD, when technically possible, in order to check if potentially large logarith- 
mic corrections to lattice artifacts predicted by perturbative analysis appear. 

4.4 Symanzik improved lattice actions 




(4.16) 




(4.17) 



If Symanzik's conjecture is true it practically follows that it is possible to find a 
Symanzik improved lattice action such that Sf s = 0, i.e. for this action there are no 
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lattice artifacts 0(a p ). The conjecture is generally accepted for AF theories, but I 
should mention that a rigorous proof of the existence of a Symanzik improved lattice 
action for any theory (including </> 4 ) to all orders PT d is, to my knowledge, not 
complete. But there is an all order proof by Keller (j 1993|) for the existence of Symanzik 
improved actions for cf>i and QED in the framework of a continuum regularization 
(using flow equations) IS A subtle point is that the continuum limit of lattice 0| 
theory is probably trivial i.e. a free theory. The renormalized coupling goes to zero 
as c/ \n(a/j,) and hence the continuum limit is actually reached only logarithmically! 
Treating the renormalized coupling g effectively as a constant for a range of cutoffs one 
has for small g a perturbative Lagrangian description of the low energy physics, and 
in this case the Symanzik effective Lagrangian describes the leading cutoff corrections 
to this. 

An important ingredient of a lattice proof (to all orders PT) would presumably 
need a proof of the small a expansion of an arbitrary ^-loop Feynman diagram of the 
form 

oo I 

F(p,a) ~a-"J2 an J2( lnam y F ™(p), (4-18) 

n=0 j-=0 

which we are quite confident holds and hence often stated in the literature, but which 
again has not, to my knowledge, been proven for £ > 2. 

Exercise 4.2 As an example of an expansion of the form (|4.18(l . consider a lattice §\ theory 
with free propagator l/[R(k, a) + m 2 ] with 

R(q/a,a) = a~ 2 r{q) , r(q) ~,=o q 2 + c(q 2 ) 2 + d^q* + ... 

Show that the tadpole integral Ji = J_J) a d 4 & a ) +m 2 ]~ 1 has an expansion of the form 

Ji ~ a~ 2 {ro + a 2 m 2 [ri + si ln(ara)] + a 4 m 4 [r2 + S2 ln(am)] + 0(a 6 )} 
with S2 = if c = d — 0. 



4.5 On-shell improved action for pure Yang-Mills theory 

With the insight gained from our previous discussion we are now prepared to consider 
improved actions for Yang-Mills theory in 4 dimensions ( Luscher, 1984[ ). As by now 
familiar, the first step is to classify the independent (up to total derivatives) gauge 
invariant operators of dimension 6, which are scalars under lattice rotations; there are 
three such operators (Wcisz, f983): 

1 = Y J tvD^D^ , 2 = J2 teD^pD^Fvp , 3 = ^ tTD„F„ p D v F vp , 

H,i/ fillip ^i^SP 

(4.19) 

2 Symanzik's published papers dealt with lower orders PT probably he considered the generalization 
straight forward. 

3 Here the improvement refers to effects involving the UV cutoff A occurring in the definition of 
the bare propagators. 
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(1) 



where D li F vp = d^F up 
the form 



(2) (3) 
Fig. 4.3 4 and 6-link closed curves on the lattice 

go [Ap, F wp \. Candidates for Symanzik improved actions have 



5i mD = y^c»(ffo) ^2 £( Ci 



C(C) 



2 



(4.20) 
(4.21) 



where U(C) is the ordered product of link variables around the closed curves C, and 
Si are sets with a given topology e.g. So = the set of plaquettes, 
Si = the set of 2 x 1 rectangles, 52 = the set of "twisted chairs" , 
5.3 = the set of "chairs" , as depicted in Fig. 14.31 

Identifying the U(x, fx) with phase factors in the continuum associated with the links 
as in (|3.22[) (with CV- replaced by goA p ), the classical small a expansion of the local 
lattice operators is given by 



ceSi,xec 



a^ZitxF^F^v + a 6 2J 



PijOj(x). 



(4.22) 
(4.23) 



We need only 4 lattice operators to represent the 4 continuum operators of dimension 
4,6 appearing in the effective action. We chose the sets of curves with smallest perime- 
ter < 6 mentioned above, but many other choices are admissible (and have appeared 
in the literature) . In order that the coefficient of F 2 in the classical expansion has the 
usual normalization the coefficients must satisfy 



co(So) + 801(50) + 8c 2 (.g ) + 16c 3 (g ) = 1 • 
Further one finds that improvement of the classical action requires 



(4.24) 



4 For this computation it is convenient to chose an axial gauge, and it is sufficient to consider 
Abelian fields. 
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Cl (0) = -1, c 2 (0) = c 3 (0) = 0. (4.25) 

On the other hand improvement of on-shell quantities only requires two conditions 
among the coefficients, e.g. for the static 2-quark potential at tree level one finds 



V(R) = ~c£L 



I : :i( Cl (0)- C2 (0)-c 3 (0)-^ R2 



(4.26) 



Spectral quantities in finite volumes, e.g. with twisted boundary conditions, can also 
be computed, and reproduces the improvement condition above and in addition the 
condition C2(0) = (Luschcr and Weisz I985&, f 985c). There are no other independent 
relations. This is as expected because on-shell the operator O3 can be dropped from 
the effective action since it vanishes when using the equations of motion. For on-shell 
improvement we conclude that we can choose 

c 3 (ffo)=0 for all g . (4.27) 

The coefficients Ci can be computed to higher orders in PT. To 1-loop they have been 
computed from the same observables as mentioned above ( |Luscher and Weisz, 1985a| ). 

In principle the Ci(go) could be computed non-perturbatively by demanding cutoff 
effects to vanish from some spectral levels e.g. requiring the J PC = 2 ++ states to be 
degenerate. This has not been done yet, the reason being that the physical goal is QCD 
and to achieve full 0(a 2 ) improvement in that theory is much more difficult because 
there are more dimension 6 operators to be taken into account, in particular those 
involving products of 4 quark fields. Never the less improved gauge actions, based on 
various considerations, are used in practical large scale numerical simulations. 
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As discussed in Hernandez' lectures, there are many ways of putting fermions on the 
lattice, each having their own particular advantages and disadvantages. In this section 
I will only discuss lattice artifacts with Wilson fermions. From the tree level coupling 
of quarks to one gluon 

~SoT a { 7A1 -ia(p+P% + 0(a 2 )} , (5.1) 

we immediately see that there are 0(a) off-shell effects, and some persist on-shell. 
The on-shell 0(a) improvement program ( |Luscher et at, 1997[ ) proceeds on the same 
lines as for the pure gauge fields in the last section. The first step is thus to classify 
the independent local gauge invariant operators of dimension 5 that can occur in the 
effective Lagrangian C\\ 

i 

Oi = goipia^F^ip , cr M „ = - [7^,7^] , 

O3 = mgltrF^F^ , 4 = imp (d^ - D^j -f^ , 5 = m 2 iptp 

On-shell we can use the equations of motion ( 7 D + m)ip = to derive relations 

Oi-O 2 + 2O 5 -0, O4 + 2O5-O, (5.5) 

which can be used to eliminate O2, O4. A Symanzik-improved action should then be 
constructible by adding a linear combination of lattice representations of 0\, O3, O5 

8S = a 5 J2 {ci(go)Oi(x) + c 3 (g )d 3 (x) + c 5 (g )O 5 (x)} , (5.6) 

X 

to the Wilson fermion action. Now 03,05 are already present in the original lattice 
action, so adding these terms merely corresponds to a rescaling of the bare coupling 
and masses by terms ~ 1 + 0(am). i.e. they can be dropped until we discuss renor- 
malization. We conclude that the on-shell improved Wilson action has only one extra 
term 

Simp = S w + a 5 '^csw{g )^{x)-<J IJ , v F IJiV (x)'il){x) , (5.7) 

X 



(5.2) 
(5.3) 
(5.4) 
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Fig. 5.1 Lattice representative of goF^ 

where the lattice representative F^ v (x) of goF flu (x) depicted in Fig. I5.1l has a "clover- 
leaf form" . This action was first written down by Sheikholeslami and Wohlert (1985). 
The coefficient csw(<?o) is known to 1-,2-loop order of PT for various gauge actions. 

Two comments are in order here. Firstly as for any improvement, 0(a) improve- 
ment is only an asymptotic concept; it could happen that at some go we have bad luck 
and 0(a 2 ) effects are bigger than 0(a) effects for some quantities. Secondly, a rather 
nice feature is that the extra overhead CPU cost in simulations of adding a SW term 
is not very substantial. 

5.1 0(a) improvement of operators 

In order to eliminate lattice artifacts in correlation functions involving composite oper- 
ators, we must also improve the operators themselves, by adding local terms of higher 
dimension having the same quantum numbers. For example for the axial iso-vector 
current A® the dimension 4 operator (A*)i in the effective Lagrangian description 
could involve terms 

O a ^ =?(A,-K)<V,75^, (5.8) 

= 9 M ?75 \r a ^ = d,P a , 3 Q M = mAl . (5.9) 

On-shell one can eliminate in favor of 0%^, O^; moreover is just a renormal- 
ization of the original operator. An ansatz for an improved lattice bare operator is 
then 

{AD^Al+acj^^ + dDP*, (5.10) 
and similarly for other operators e.g. 

pa = pa ^ { ya^ = ya + ^ ^ + ^) ^ ^ _ (5 U) 

The improvement coefficients cv,ca appearing above are O(g^). 

5.2 Mass independent renormalization scheme 

In lattice QCD it is often advantageous to use a mass independent renormalization 
scheme as discussed in sect. 2. Without 0(a) improvement this would involve renor- 
malization constants independent of the quark masses and take the form 
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Sr = So%(SoW) , (5.12) 
m R = m q Z m (gl,afx) , m q = m - m c . (5.13) 

To obtain correlation functions at physical distances to approach their continuum 
limit at a rate 0(a 2 ) with improved Wilson fermions, we must modify the form of 
the renormalized parameters (and thereby account for the terms 03,05 we dropped 
previously) in order to avoid uncanceled 0(am q ) effects ( Luscher et al., 1997 1. We do 
this by introducing modified bare parameters 

So = So (! + bgam q ) , (5.14) 
m q = m q (1 + b m am q ) , (5.15) 

and define renormalized parameters through 

Sr. = 5o^s(5oi a M)) ( 5 - 16 ) 
m R = rh q Z m (g%, afi) . (5-17) 

We can compute the new improvement coefficients in perturbation theory, e.g. from 
the "pole mass" at tree level 

mp = — ln(l + amo) — m q — —am q + . . . , (5.18) 

we deduce 

bm = -l + O(g 2 ). (5.19) 

The lowest contribution to the coefficient b g is 0(<?o) and first obtained by Sint and 
Sommer (1996) from the 1-loop computation of the SF coupling with fermions at fixed 
z = m q L: 

9sf(L) = So + So [ 2& o ML /a) + C(z) + Nfkam q + ...], (5.20) 
from which one deduces b g — Nfkg^ + . . . . 

Similar factors are required for renormalized composite operators 

= Z A (1 + b A am q ) (A a j)^ , (5.21) 
P£ = Z P (1 + b P am q ) P a . (5.22) 

The coefficients 6a, bp = 1 + 0(<?q), which were first computed to 1-loop order by 
( jGabrielh et al, 1991[ ), don't depend on renormalization conditions - the latter are 
applied at m q = 0. 

5.3 Determination of RGI masses and running masses 

We recall from sect. 2 that renormalization group invariant masses are scheme inde- 
pendent. Running quark masses can be defined non-perturbatively in various ways; a 
most efficient way makes use of the PCAC relation: 

Z A (dA a O) 

m(M) = 2Z p (a0(P°O) ' (5 ' 23) 

where we have not yet specified the source O (the running mass should be practically 
independent of this) nor the precise definition of the expectation value. The running 
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of a PCAC mass is determined by the running of the renormalization constant of the 
pseudoscalar density Zp. 

One scheme for computing scale dependent renormalization constants e.g. Zp, is 
the regularization independent (RI) MOM scheme described for couplings in sect. 3.1. 
Originally introduced by ( |Martinelli et al., 1994| ), it is now quite popular and used by 
many collaborations. Progress using this scheme was reported by Y. Aoki at LATT09 
( |Aoki, 2008] ), flSturm et al., 2009] ) and will be covered in these lectures by Vladikas. 

Again an alternative is to apply hnitc size recursion techniques here too, and defin- 
ing Zp running with the volume. Then in the continuum we can define an associated 
step scaling function ap through 

m(l/L) Zp{2L) _ 2 

4M = ^y = ap{uh u = g {L) - (5 - 24) 

In perturbation theory we have 

fip(u) = 1 -d ln(2)u + O(u 2 ). (5.25) 

After having chosen a definition of Zp on the lattice for a given go,L/a, we can 
proceed to compute op{u) via 

<7p(w) = lim Ep(u, a/L) , (5.26) 

a— >0 

£ P (u, a/L) = — — rj -^- g2 (i)= „ . (5.27) 

Zp{go,L/a) 

Suppose we have already determined the step scaling function for the running coupling 
(as described in sect. 3) at say 8 points: 

L k = 2- k L max , u k =f{L k ), fc=l,...,8. (5.28) 

Then at each step k we measure T,p(u k ,a/L) for a sequence of values of L/a and 
extrapolate these to the continuum limit: 

Ep(u fe ,a/L) = crp(u fe ) + 0(a p ) , (5.29) 

where p — 2 for the O(a) improved theory (up to boundary terms). Thereby we have 

M M 



J\op{u k ). (5.30) 



m(l/L 

max 

) m(l/L 8 ) 

Since L% is a very small physical length we can safely use PT to determine the first 
factor _ jM-f , . To relate the RGI mass M to a low energy scale it then remains to 

m(l/L s ) ° J 

compute I max and m(l/L max ) in physical units e.g. f v . 

5.4 The Alpha Collaboration project 

It has been an Alpha Collaboration goal since ~ 1998 to measure running couplings, 
running quark masses and renormalization constants in QCD to a high precision. The 
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program uses 0(a) improved Wilson fcrmions in the SF framework. As we have tried to 
emphasize this procedure is comparatively clean and care is taken to control systematic 
errors at each stage, but it involves a lot of preparatory work and attention to details. 

Firstly there are extra improvement coefficients needed to cancel 0(a) boundary 
effects, even for the pure gauge theory ~ J gA F^ k , f gA F i | . This involves adjusting 
weights for the plaquettes at the boundary 

»(P) = I I 2 J Cs(ff ° ),f ° r ?f A v A ). ( 5 - 31 ) 
t c H<?o) i° r P temporal touching A J 

with c Si t = 1 + 0(oq). There are similar coefficients c s ,c~t for fermionic boundary 
terms. Moreover there is a ^-coefficient for the renormalized boundary operator £r = 
Z^(g^, afi)(l + b^arriq) although these factors usually cancel in ratios defining observ- 
ables of interest. It has been observed that these coefficients have small coefficients, 
and thus in practice it is considered safe to use the 1-loop perturbative approximation 
instead of determining them non-perturbatively. 

Many of the steps must be carried out in a definite order. For example the SF 
running coupling is defined as before and at zero quark masses m q = 0. This means 
that before measuring the coupling the critical mass m c must be determined but this 
in turn depends on c CS w- In this project m c is defined through the vanishing of the 
PCAC mass - other definitions e.g. by vanishing pion mass will differ by 0(a 2 ) in the 
improved theory. 

In an improved theory we expect ( |Luscher et al, 1997[ ) 

((A R )«{x)J a ) = 2m K {P&(x)J a ) + 0(a 2 ) (5.32) 

for arbitrary sources J a . In the framework of the SF we can take J a = O a as defined 
in (I3.3ip . In terms of the bare correlation functions 

f A (xo) = -\{Al{x)O a ) , (5.33) 

fp(x ) = -±(P a (x)O a ), (5.34) 

we then define improved bare PCAC masses by 

i a n n'\ + d*)f Al {x Q ) 
m(x ,6 k ,C,C ) = —— — . (5.35) 

Finally renormalized SF PCAC masses are given by 

- ,r\ (l + b A am q )Z A (g ) 

m SF (L) = —— -rj-rm, (5.36) 

(1 + bpam q )Zp(g ,L/a) 

where to complete the definitions we must specify the arguments of m. 

It is clear that the coefficients ca, c csw can be determined by demanding that m 
is independent of the sources up to terms 0(a 2 ). This can be done by appropriately 
choosing independent configurations of xo, Ok, C, C . One can proceed as follows: noting 
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that m(xo) is linear in ca, m(xo) = r(xo) + cas(xq) we can form from two choices of 
boundary conditions (Ci, C[), (C 2 , C' 2 ) at say 9k — 0, a linear combination 

as 1 \ 1 \ 1 \ I TOl (yo) - m 2 {yo)] q ^ 

Mi(x ,y ) = mi i - «i(zo) 7—^ t^-t^ (5.37) 

si(y ) - 52(2/0) 

which is independent of ca and equal to mi(xo) up to 0(a 2 ) corrections (in the im- 
proved theory). We can now give a condition determining c csw e.g. 

= M 1 (T/4,3T/4)-M 2 (T/4,3T/4) (for M^T/2, T/4) = 0) . (5.38) 

A condition such as the latter in brackets is necessary to complete the specification; to 
ensure that the results are relatively insensitive to the precise choice one should check 
that the dependence on m is weak for m ~ m c (which is only known approximately at 
this stage). Once c csw has been determined, we can specify a computation of ca e.g. 
by comparing values of m at C = C = and varying Ok- 

The data for c csw ,ca thus obtained are usually fitted in the measured range to 
(rational) functions of go which incorporate the known perturbative coefficients. These 
functions can now be considered as definitions of the improved theory, which can then 
be used in other simulations (not necessarily in the SF framework). 

It is important to appreciate that c csw ,ca have 0(a) ambiguities; there is no way 
to remove these outside of PT. Once these coefficients have been specified, we can use 
the vanishing of mpcAC for a determination of m c . Again the precise value depends 
on all specific definitions involved. 

Knowing c csw , m c we can now compute g~s F (L) and obtain the step scaling function 
in the continuum limit. 

Next one can compute the pseudoscalar density scaling function ap as described 
in sect. 5.3 in a given SF scheme e.g. 

Zp(go,L/a) = j- p ^y2) \T=L, m =m c ,c=c>=o,e k =i/2 , (5.39) 
h - -^(O la O a ) , (5.40) 

and c is chosen st Zp\ go= o = 1. 

Finally one can compute the vector isovector current renormalization constant Zy 
using the WI (the derivation is left as an exercise) 

(1 + byamJZvMxo) = fi + 0(a 2 ) , (5.41) 

where 

/vK) = ^((^)S(^t), (5-42) 

0^ t = e abc O b O' c , (5.43) 



and for Za- 
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(1 + b A am q )Z A f AlAl (x ,y ) = fi + 0(a 2 ) . 



where 



fAA(x ,yo) 



L 6 



Y,e abc (A a (x)A b (y)O: 



(5.44) 



(5.45) 



the conditions for determining Z A , Zy being at m q = 0. No explicit results are shown 
here but the program outlined above is practically complete for Nt = 2 species of dy- 
namical fermions ([Delia Morte et al. 2005al [l gll2005&|), (Delia Morte et al, 2005 c| , 
(Delia Morte et al, 2005^ ), and for N { = 3 flAoki et al, 2009[ ). In particular the SF 



coupling is measured over a wide range of energies ( Delia Morte et al., 2005 a| with a 
result similar in quality to that for N{ = 0. Finally some results for the Alpha program 
are recently available for Nf = 4 flTekin et al., 2010] ). 



6 

Other improved actions 



Many types of improved actions have been used in large scale simulations. As men- 
tioned before this is also useful for numerically verifying universality. Often due to the 
pressure of completing the simulation routines before a new supercomputer is deliv- 
ered, pragmatic choices of the action to be used have to be made. For example the 
CPPACS collaboration decided to employ the Iwasaki action ( jlwasaki, 1983[ ) which 
has the form (|4.20[) with cq = 3.648, c\ = — 0.331, C2 = c 3 = independent of go, 
but optimized for values of go in the range to be simulated (e.g. by demanding good 
rotational symmetry properties of the static potential). 

6.1 Perfect actions 

An action used by Hasenfratz, Hasenfratz and Niedermayer (12005| is based on Wilson's 
renormalization group approach (see Niedermayer's review (|1998|) ). The intriguing 
realization is that in the huge class of lattice actions which have the same universal 
continuum limit there are "perfect actions" for which the physical quantities that can 
be measured have no lattice artifacts at all even when the correlation lengths are 0(a)! 

To understand how this comes about consider an RG transformation in configura- 
tion space of a pure gauge theory of the form 

e -P'A'on = J [du] e -p[A(u) + nu,v)] (fU) 

with kernel 

T(U,V) = -^J2 mi [V(n B ,IJ)QHn B ,fj)-N^ , (6.2) 

where Q(ub,Ij) is a sum over products of link variables U on paths from 2ns to 
2(ns+A) on the original lattice (e.g. a sum over staples). For physical quantities which 
can be measured with action /3'A' we will get the same results as those measured with 
ft A. The main problem however is that we are now dealing with an infinite coupling 
parameter space. 

The situation is simplified for the case of asymptotically free theories with one 
relevant direction labeled by ft, with the critical surface where the correlation length 
diverges given by ft = oo. If we start the RG transformation with ft very large then 
we expect that A' will have an expansion of the form 



A'(V) = A (V) + -A 1 (V) + 



(6.3) 
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C 2 




Fig. 6.1 RG flow 



(see Fig. 16.11) where 

A Q (V) = min [A(U) + T(U, V)} . (6.4) 

Performing this minimization with A (U) on the critical surface /3 = oo we stay on that 
surface. Repeating it infinitely many times we finally reach the "fixed point action" 
A FP associated with the kernel T. It satisfies the equation 

A FP (V) = min [A FP (U) + T(U, V)] . (6.5) 

This action is classically perfect; there exist solutions of the classical field equations 
(e.g. instanton configurations) with no lattice artifacts. If one starts the RG transfor- 
mation close to the FP one obtains an RG trajectory along which one has very small 
lattice artifacts even for very small correlation lengths. 

The proposal of Hasenfratz and Niedermayer (|1994[) is to simulate the action (3A FP 
for f3 large. In practice it is too time consuming to determine A FP configuration by 
configuration. Instead one does it for a large number of configurations and parameter- 
izes the results with a relatively small number of loops. Also they tune the parameter 
K in (|6.2[> to weaken the the strength of the spread of loops in lattice units. 

Fermions can be included along similar lines. A very nice feature is that fermion ac- 
tions which are obtained by blocking from a continuum action which has chiral symme- 
try have Dirac operators which satisfy the Ginsparg- Wilson relation ( Ginsparg and Wilson, 1982| 

{V^ 5 } = (l + s)V l5 V, (6.6) 

(or some mild local generalization thereof), where s is a parameter |s| < 1. Consider for 
example lattice fermions x obtained from blocking a free massless continuum field tp. 
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The lattice Dirac operator can be obtained by a minimization procedure over classical 
fields: 

Xn^nn'Xw = min {ipDip + (x - i>^) (x ~ , (6-7) 

where 

D xx> = h^) xx , , (6.8) 

and 

w nx = 1 if x G block n 

= otherwise. (6-9) 

The minimizing field is given by 

Mx)=A~ 1 ^x, Mx)=XloA-\ (6.10) 

with 

A = D + U3 l u. (6.11) 



Exercise 6.1 Show that the lattice Dirac operator 

V = 1 -cjtTV (6.12) 
satisfies the GW relation (16.61) with s = 1. 



Moreover Hasenfratz et al show in a very illuminating paper (12006j) the general 
result that lattice actions induced by a RG procedure inherit all the symmetries of 
the continuum theory, and they give a general procedure which delivers the corre- 
sponding symmetry transformation on the lattice e.g. for the U(l) axial continuum 
transformation 

5ip = iejsip , Sip = ietpj5 , (6.13) 
one obtains Luscher's lattice transformation ( jLuscher, 1998[ ) 

S x = in 5 {l-V) x , Sx = tex(l - £>)75 . (6.14) 
6.2 Neuberger's action 

There are also classes of Dirac operators satisfying the GW relation which are not ob- 
tained by RG considerations e.g. Neuberger's overlap massless Dirac operator ( |Neuberger 1 998a, 
I19986|> 



D - (1 + S) 
D N — 



l-A(A^A) 1/2 , (6.15) 
a L J 

A = l + s- aD w , (6.16) 



which is discussed in Kaplan's lectures (see also Niedermayer's review (|1999j) ). Sim- 
ulations with associated actions are highly desirable because of their excellent chiral 
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properties, however at present they are very (CPU) expensive and hence only pur- 
sued by a minority of groups such as JLQCD QOhki et al., 2009[ ). Here I only want to 
summarize a few points concerning their renormalization and their lattice artifacts. 

The Feynman rules are straightforward to derive but a bit more complicated than 
other actions because the action is not ultra-local. The free massless propagator is 
rather simple and takes the form (p^ — a -1 sm(ap^)) 

So{p) = 2(TT^) {-^^ + b ^ + ^ 

b(p) = il±^ - , u (p) = ■ (6.18) 

a 2 

The quark-antiquark 1-gluon vertex V(p, q) has a relatively familiar structure, however 
the quark-antiquark 2-gluon vertex involves also an integral over the product of two 
such vertices ~ J V(r, p)V(r, q)K(p, q, r) with some kernel K . With such a structure it 
is not immediately clear that the Reisz conditions for power-counting are met, however 
Reisz and Rothe (2000) have proven the renormalizability of Neuberger fermions with 
2 species of massless quarks. 

As explained in Kaplan's talk the massless GW action can be decomposed in a 
sum of left and right-handed parts 

^V^j = i> L V^ L + i> R Vi]} R , (6.19) 



with 

ipR/L =P±tp, 
where the projectors are given by 

P± = 5(1 ±7s) , 
75 = 75(1 - aD) . 



^r/l = , (6-20) 

?± = \{1±%), (6.21) 

(6.22) 



It follows that scalar and pseudo-scalar densities are naturally defined by 

S = ^ L ^ R +^ R i, L =^(l-^D)i>, (6.23) 

P = Jj l ^r - ~ip R ipL = ^75 (l - ~X>) V , (6-24) 

which transform into one another under the chiral transformation 

8ipL = -i>L, 6ip L =ip Ll 5ipR = = 5ip R . (6.25) 

Mass terms are naturally introduced into GW actions via 

V -> V + mS . (6.26) 

The mass is renormalized multiplicatively as in the continuum. The pattern of renor- 
malization of quark bilinears and 4-quark operators is also the same as in the contin- 
uum, and for Neuberger fermions computed to 1-loop perturbation theory by Capitani 
and Giusti flUO). 

1 their proof uses the fact that for Nf > 2 GW lattice fermions possess an exact chiral (flavor 
mixing) symmetry. 
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Exactly conserved flavor vector currents which do not require renormalization have 
been constructed |Chandrasckharan, 1999), (Kiku kawa and Yamada, 1999[ ), (Hascnfra tz et al., 2002). 
Since these are relatively complicated other currents of the form t/^ 5 ( 1 — ip are 
usually used as observables. They also have nice chiral transformation properties but 
are not exactly conserved and hence require finite renormalization (Alcxandro u et al, 2000[ ). 

The discussions of lattice artifacts are again based on accepting the validity of 
arguments based on Symanzik's effective action. Firstly for to = there are no 0(a) 
artifacts for spectral quantities since in this case any such effects would be described 
by operators of dimension 5 having symmetries of the lattice action. No such operator 
can be constructed for GW fermions (a SW-like term breaks the chiral symmetry). 
When to 7^ the GW action is invariant under the transformation tp — > 757/', tp — > 
—7/175 together with m — > —to. This symmetry is however not quite sufficient to 
disqualify O(am) artifacts in correlation functions since these transformations change 
the measure for topological non-trivial configurations. One can argue however that in 
large volumes one can, for spectral quantities, restrict attention to the topologically 
trivial sector and hence also expect no O(am) artifacts. 

Matrix elements of composite operators mentioned above do however have 0(a) ef- 
fects. The improvement of bilinears are obtained by the prescription (to — 0) (Capitani et al. 1999, 
12000)) : 

(^rV0i = ^(l-|£)r(l-|£>)v,. (6.27) 

As a last point I would like to mention that Schrodinger functionals for overlap 
fermions have been constructed by Luscher ( |2006j) and in a different way by Sint (for 
Nf even) (|20076|) . In order to obtain the natural SF boundary conditions (|3.28p in the 
continuum limit, Luscher simply modifies the Neuberger operator so that it satisfies 
the GW relation up to terms on the boundary 

{ l5 ,V} = (l + s)aV l5 V + A B (6.28) 

where AbiI>{x) — , unless xq = 0, T. A simple proposal is 



2> = £±i> 



1-^(^ + 75^75) 



(6.29) 



with 



U = A (A' A + caP) 1/2 , (6.30) 

Plj){x) = - {S Xo a P-1p(x)\ Xo = a + $x .T-aP+lp{x)\x =T-a} , (6.31) 

a 

where c is a parameter which can be tuned for O(a) improvement, and has been 
computed to 1-loop order PT by Takeda (|2008"1) . 



6.3 Twisted mass lattice QCD 

As the last topic I would like to mention Wilson twisted mass lattice QCD; and refer 
the reader to excellent reviews by Sint (12007a|) and Shindler (|2008|) . For the case of 
two flavors (N[ = 2) the action takes the form 
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Stm = J dxdx X d tmX , (6-32) 
Dtm = D w + mo + i^ q 75T 3 • (6.33) 

It appeared first in papers by Aoki and Gosch (1989, 1990) who were considering 
the extra twisted mass term as an external probe to study the phase diagram of 
Wilson fermions. In these works it was realized that the associated hermitian operator 
Q = j^Dtm has no zero modes for /x q > 0: 

Q = Qw + i/V" 3 ! an d Qw — Qw Q^Q = QwQw + A*q ■ (6.34) 

Thus adding a twisted mass term provides a local field theoretic solution to the problem 
of treating exceptional configurations in quenched simulations. The TM formulation 
also has some advantages in dynamical HMC simulations. 

Soon it was realized (|Frezzotti et al. 2001 al IFrezzotti et al. 2001 b\ that TMQCD 
is a viable alternative regularization of QCD, and it is, for an additional reason that 
I will soon explain, now used in some large scale simulations, This may at first sight 
seem strange because parity is apparently strongly violated, but this is not the case. 
To see this consider first TMQCD in the formal continuum limit 

Stmqcd = J X (jD + TO q + i/i q 7 5 r 3 ) x ■ (6-35) 

Now making a change of variables in the functional integral 

1> = V{u)x, ^ = W(u), (6.36) 
V(ui) = exp ( iujjs — ) , with tanw = —3- . (6.37) 

V 2 ) TO q 

the action takes the usual form 

Stmqcd = J^(r/D + M)i/> , (6.38) 

with 



M = ^jml + M 2 . (6.39) 

Hence continuum twisted mass QCD is formally equivalent to QCD. The symmetries 
in one basis are in one-to-one correspondence with symmetries in the other basis. The 
■0's are called the "physical basis" because in this basis the physical interpretation is 
clearer and x's are called the TM basis. 

The formal equivalence can be made more rigorous dFrezzotti et al., 2001 q| ) by 
considering a lattice regularization of TMQCD which preserves a chiral symmetry such 



2 These are configurations which have exceptionally small eigenvalues of Qyy and hence give huge 
contributions to some correlation functions e.g. (P a (x)P a (y)) in quenched simulations. 
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as GW fermions. In such a case the regularized correlation functions in the different 
bases are equal: 

mM(M,0) = (O[x,x})(m cl , IM )- (6.40) 

Operators in one basis are equivalent to associated operators in the other; this requires 
some familiarization e.g. for fermion bilinears in the TM basis define 

K = Xl» yX , similarly A«, S a , P a and S° = XX, (6.41) 
and in the physical basis 

K = V^y^ - similarly A%, S\V a and 5° = ^ . (6.42) 
Then these are related by 



A* = cos ujA^ + e 3ab sin u , a =1,2, (6.43) 
= A», a = 3, (6.44) 



and similarly for V£, and 



V a = P a , a =1,2 (6.45) 
= coscjP 3 + 2ismujS° , a = 3, (6.46) 
5° =cosu;5' + 2isinwP 3 . (6.47) 



The familiar PCVC and PCAC relations 



W = 0, (6.48) 
d„A% = 2MT a , (6.49) 



translate in the TM basis to 



8^ = -2^e 3ab P b , (6.50) 
d^A* = 2m q P a + ifi q S 3a S° . (6.51) 

For the case of a regularization which breaks chiral symmetry, such as Wilson TM 
fermions, one expects that near the continuum limit 

(0R[X,X] W q ) = {Ok[^W°m% + 0(a) . (6.52) 

The lattice artifacts depend on m q and on the twist angle to. In this sense ordinary 
Wilson fermions are just a special case of the WTM regularization of QCD. 

To discuss the lattice artifacts the determination of the structure of the Symanzik 
effective action again plays a central role. This first requires listing the symmetries of 
the WTM action; among these are 

1. gauge invariance, lattice rotations and translations, and charge conjugation 
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2. P^-symmetry which is conventional parity V combined with a flavor rotation: 
V% : x(x) -► ir a P X (x) , a = 1,2 7>x(x ,x) = 7ox(z , -x) 

3. and "P = 7? x {^ q -> -//J 

4. similarly for time reversal 

5. also the Uy(l) subgroup of SUy(2) is maintained 

With this knowledge one can again show that the only terms in the on-shell effective 
action (apart from ones that can be incorporated in the renormalizations of the bare 
parameters) is the SW term 



(6.53) 



The 0(a) improvement program of the Wilson fermions can be extended to the case of 
arbitrary twist angle uj =^ . In general more o-coeflicients are required e.g. for the mass 
independent renormalization one has to introduce improved bare mass parameters 



rn. 



q — m q (l + b m am q ) + l m a^\ , 
Mq = Mq(l + bfj.amq) , 



(6.54) 
(6.55) 



and for the rcnormalized axial current 



(A K )l = Z A {1 + b A am q )[A« + ac A d^P a + a f i q b A e 3ab V. b 



(6.56) 



An obvious simplification of the regularization at maximal twist angle u) = tt/2 
was realized early on, but TMQCD only became more popular as a practical regular- 
ization after the realization by Frezzotti and Rossi (|2004[) of a related property called 
"automatic 0(a) improvement". This is the property that properly renormalized phys- 
ical expectation values have only 0(a 2 ) artifacts provided that the untwisted mass is 
tuned to its critical value, e.g. by tuning the PCAC mass oc (d^A^O) to 0. The original 
demonstration ( |Frezzotti and Rossi, 2004] ) of this was not complete, but the discov- 
erers and many other authors developed independent proo fs (|Sharpe and Wu, 2 005), 
( Shindler, 2006] ), ( jSint, 20060 , (jAoki and Bar 20041 120U6)) . ( [Frezzotti et al, 20060 . 

Working in the TM basis with 



So = J d 4 xx [lD + ifi q -/ 5 T 3 ] x 



we have 



(Or) = (0 R ) cont - a(0 R Sf) cont + a(Of) cont + 0(a 2 ) 
One proof is based on the following symmetry of Sq: 



x ■ 



(6.57) 
(6.58) 

(6.59) 

(6.60) 



IZl alone is not a symmetry of WTM, but Hi, x D is, where V essentially measures 
the parity of the dimensions of operators 
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V : U{x,n) -+ U\-x - afi, n) , (6.61) 
x(x ) e 3 ^ 2 X (-x) , x(x) -> e 3 -/ 2 x(-x) . (6.62) 

Consider now operators C± which are even(odd) under 72.5. Now S 1 ^ is odd under 1Z\ ; 
also (using I?) the effective operators are odd(even) under IZl, and so it follows 

(O+r) = (0 + n) cont + 0(a 2 ) , (6.63) 
(O-r) = 0(a) . (6.64) 

One can extend the TM regularization to include non-degenerate quarks and also 
further flavors in various ways. 

Another advantage of using the TM formulation is in the simplification of some 
renormalization properties of composite operators. For example in the computation of 
f n with Wilson fermions from the correlator {•A^a{x)V^[y)) one needs but using 
the equivalence 

« (x)^(y))(M R ,o) =s (V 2 (x)Pi(y)} iMK ,,/2) (6.65) 

no renormalization constant for the current on the rhs is needed. 

Another example is the 4-fermion operator appearing in the description of AS" = 2 
transitions. In the physical basis 

°tv-A){V-A) = E " ^)rf) 2 (6-66) 

is a sum of two parts, one parity even Ovv+^^l an d the other parity odd —Ova+av- 
As mentioned in sect. 2.7 for Wilson fermions the parity even operator mixes with 
many other operators. However with WTM regularization at maximal twist in the u, d 
sector and Wilson for the s-quarks one obtains 

Ovv+aa - V a+av ■ (6.67) 

Using C, P symmetry one can show that Ova+av renormalizes diagonally, which 
simplifies the problem considerably QPena et al, 2004[ ). 

In conclusion, TMQCD is a competitive practical regularization of QCD, especially 
at maximal twist. The price to be paid however is a breaking of parity and flavor 
symmetries. This is usually at the 0(a 2 ) level but these effects can in some situations 
be uncomfortably large. 
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